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SECTION  1 


INTRODUCTION 


Non-reflecting  boundaries  for  transmitting  transient  wave  energy  in  an  unbounded  elastic  solid  have 
many  applications,  including  ground  shock  analysis,  earthquake  analysis,  and  quantitative 
non-destructive  evaluation.  For  example,  in  a  dynamic  analysis  of  a  fully  or  partially  buried  structure, 
the  structure  is  of  primary  concern,  but  a  certain  amount  of  the  surrounding  soil  is  typically  included 
in  the  analysis  to  treat  the  soil-structure  interaction.  Unfortunately,  spurious  reflections  intrude  into  the 
calculations  when  waves  scattered  from  the  structure  reach  a  fixed  or  free  soil  boundary.  This  forces  the 
analyst  to  keep  the  boundary  far  from  the  structure,  thereby  requiring  the  inclusion  of  an  large  amount 
of  soil.  Thus,  non-reflecting  boundaries,  which  may  be  placed  relatively  close  to  or  directly  on  the 
embedded  structure,  are  needed. 

Because  of  the  broad  spectrum  of  frequencies  that  need  to  be  included  in  many  transient  dynamics  anal¬ 
yses,  it  is  not  efficient  to  employ  Fourier  transform  (time-harmonic)  methods.  More  importantly,  when 
the  near-field  medium  and/or  structure  respond  nonlinearly,  time-harmonic  methods  cannot  be  used; 
only  transient  methods  can  be  used.  Doubly  Asymptotic  Approximations  (DAA’s)  are  time-domain 
non-reflecting  boundary  relations  that  match  asymptotically  the  exact  elastodynamic  integral  represen¬ 
tation  for  infinite  and  semi-infinite  domains  at  both  early  and  late  times.  For  broad-spectrum  problems, 
doubly  asymptotic  approximations  are  more  robust  than  singly  asymptotic  approximations,  which  only 
match  the  exact  representations  at  one  limit. 

Currently,  most  non-reflecting  boundaries  implemented  in  transient  finite-element  analysis  programs 
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are  based  on  the  plane-wave  (early-time)  approximation(see,  e.g.,  (Whirley,  1991)),  which  is  singly 
asymptotic.  If  the  problem  of  interest  involves  late-time,  low-frequency  components,  the  plane-wave 
boundary  will  produce  incorrect  results.  However,  the  second-order  DAA  derived  during  this  study 
will  reproduce  with  high  fidelity  both  the  early-time  and  late-time  responses,  and  will  reproduce  with 
good  accuracy  the  intermediate-time  response. 


Underwood  and  Geers  (Underwood  and  Geers,  1981)  heuristically  formulated  and  Mathews  and 
Geers  (Mathews  and  Geers,  1987)  developed  further  the  first-order  DAA  for  an  infinite  domain  (DAAf" ). 
The  goal  of  the  present  study  has  been  to  formulate  and  evaluate  doubly  asymptotic  approximations 
for  transient  elastodynamic  analysis,  especially  second-order  approximations,  that  have  not  been  pre¬ 
viously  formulated.  The  specific  objectives  under  this  goal  are: 

1 .  Develop  the  first-order  doubly  asymptotic  approximation  for  a  homogenous,  semi-infinite,  elastic 
domain  (a  half-space):  DAAj^.  To  achieve  this,  a  first-order  late-time  approximation  is  derived 
and  a  systematic  operator-matching  method  is  employed  to  combine  the  early-time  and  late-time 
approximations  into  DAA(^. 

2.  Develop  the  second-order  doubly  asymptotic  approximation  for  a  homogenous,  infinite,  elastic 
domain  (a  whole-space):  DAA^.  To  achieve  this,  second-order  approximations  for  early-time 
(high  frequency)  and  late-time  (low  frequency)  response  in  infinite  domains  are  formulated.  In 
addition,  a  second-order  operator  matching  method  is  employed  to  combine  the  early-time  and 
late-time  approximations  into  DAA^. 

3.  Develop  the  second-order  doubly  asymptotic  approximation  for  a  homogenous,  semi-infinite,  elas¬ 
tic  domain  (a  half-space):  DAA^.  The  principal  challenge  in  this  effort  is  determining  a  second- 
order  late-time  (low  frequency)  approximation  for  a  semi-infinite  elastic  domain.  The  same 
operator-matching  method  used  in  the  previous  component  is  used  to  contruct  the  DAA^. 

4.  Obtain  exact  and  DAA  relations  for  the  dilatational,  rotational,  and  translational  motions  of  a 
spherical  boundary  in  an  infinite  elastic  medium  subjected  to  transient  internal  tractions. 

5.  Compare  numerical  DAAj'^,  DAA(^,  DAA^^,  and  DAA^  results  for  selected  evaluation  prob¬ 
lems  with  those  of  other  investigators.  This  involves  the  development  of  general-purpose  com¬ 
puter  programs  based  on  the  boundary  element  method.  In  addition,  comparisons  are  made  be¬ 
tween  the  analytical  modal  results  of  the  previous  objective  and  corresponding  results  generated 
with  the  general-purpose  DAA  programs. 
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How  these  objectives  are  reached  is  documented  in  the  following  sections.  The  literature  in  the  field  is 
reviewed  in  Section  2.  In  Section  3,  the  fundamental  equations  used  in  this  study  are  presented.  Sec¬ 
tion  4  addresses  the  first-  and  second-order  eariy-time  (high-frequency)  approximations,  and  Section  5 
addresses  the  first-and  second-order  late-time  (low-frequency)  approximations.  In  Section  6,  the  first- 
and  the  second-order  doubly  asymptotic  approximations  are  formulated  by  operator  matching. 

Section  7  presents  exact  and  DAA  relations  for  breathing  and  translational  motions  of  a  spherical  bound¬ 
ary  in  an  infinite  elastic  medium.  In  Section  8,  various  DAA  results  are  compared  with  results  obtained 
by  other  means.  Finally,  Section  9  states  the  conclusions  reached  in  this  study. 
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SECTION  2 


REVIEW  OF  THE  LITERATURE 


In  this  chapter  we  review  previously  published  work  in  the  field  of  transient  elastodynamics  as  applied 
to  non-reflecting  boundaries.  We  will  concentrate  on  the  theoretical  work  more  than  on  the  numerical 
work.  The  numerical  aspects,  especially  for  non-reflecting  boundaries,  may  be  found  in  Manolis  and 
Beskos  (Manolis  and  Beskos,  1988). 

2.1  ELASTODYNAMIC  THEORY. 

The  foundations  of  dynamic  elasticity  were  laid  by  Poisson,  Stokes,  Kirchhoff,  and  Love  (Poisson,  1 820), 
(Stokes,  1849),  (Kirchhoff,  1883),  (Love,  1904b),  (Love,  1904a),  (Love,  1944).  Poisson  formulated 
the  initial  value  problem  for  a  dynamic  elastic  domain  in  1819.  Poisson’s  formula  results  from  trans¬ 
forming  the  wave  equation,  cf)  =  into  an  integral  formula  that  equates  the  value  of  ^  at  a  point 

to  an  integral  taken  over  a  sphere  surrounding  the  point.  Kirchhoff  generalized  this  formula  to  equate 
the  value  of  ^  at  a  point  to  an  integral  taken  over  an  arbitrary  surrounding  surface.  Stokes  established 
the  fundamental  singular  solution  of  the  equations  of  elastodynamics  in  an  infinite  region  due  to  a  time- 
dependent  concentrated  body  force.  Love  added  to  this  foundation,  proving  that  certain  initial  velocity 
conditions  involving  continuity  had  to  be  met  for  the  integral  equations  of  Poisson  and  Kirchhoff  to  be 
applied  properly  (Love,  1904b),  (Love,  1904a),  (Love,  1944). 

Doyle  (Doyle,  1966),  using  results  from  Sternberg  and  Eubanks  (Stemburg  and  Eubanks,  1957),  deter¬ 
mined  that  the  fundamental  singular  integral  equations  could  be  simplified  by  employing  the  Laplace 
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transform.  In  addition,  Doyle  showed  that  the  initial-  and  boundary-value  problems  could  also  be  for¬ 
mulated  using  the  Laplace  transform.  Cruse  and  Rizzo  (Cruse  and  Rizzo,  1968)  extended  Doyle’s  re¬ 
sults  to  the  solution  of  general  problems  by  using  boundary  element  methods  in  the  Laplace  transform 
domain.  The  drawback  of  this  approach  is  the  difficulty  of  transform  inversion. 

If  the  boundary  integral  equations  are  not  Laplace  transformed,  they  contain  convolutions  in  time.  The 
solution  of  geometrically  complex  problems  through  the  use  of  time  convolution  is  very  computer- 
resource  intensive.  Therefore,  neither  Laplace-transformed  nor  time-convolved  boundary  integral  equa¬ 
tions  are  suitable  for  use  as  non-reflecting  boundaries  in  practical  transient  elastodynamics 
analysis  (Manolis,  1983). 

2.2  NON-REFLECTING  BOUNDARIES. 

The  objective  of  a  non-reflecting  boundary  in  scattering  analysis  is  to  completely  absorb  the  scattered 
waves  in  order  to  simulate  an  infinite  domain.  Some  early  work  on  scattering  in  an  elastic  solid  was 
done  by  Ying  and  Truell  (Ying  and  Truell,  1956),  who  formulated  a  solution  for  the  problem  of  a  plane 
longitudinal  wave  scattered  by  a  spherical  obstacle  in  an  isotropically  elastic  solid.  This  formulation, 
like  so  many  others  (Pow  and  Mow,  1972),  (Uberall  and  others,  1990),  assumes  the  incident  and  scat¬ 
tered  waves  are  time  harmonic,  which,  as  mentioned  above,  is  not  applicable  to  problems  involving 
non-linear  transient  response. 

A  second  method  used  in  scattering  problems  is  ray  theory  (Achenbach  and  others,  1982).  Ray  theory 
has  been  widely  used  in  electromagnetics,  where  it  is  often  called  geometrical  optics.  Ray  theory  em¬ 
ploys  asymptotic  approximations  that  are  accurate  for  high-frequency 
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response  (Achenbach  and  others,  1982),  (Freidlander,  1958).  Ray  theory  was  extended  by  H.  B.  Keller 
and  by  J.  B.  Keller  and  his  coworkers  from  electromagnetics  and  acoustics  to 
elastodynamics  (Keller,  1964),  (Keller,  1958),  (Ahluwalia  er  a/.,  1969). 

Recently,  Givoli  and  J.B.  Keller  have  developed  non-reflecting  boundaries  for  transient  elastodynam¬ 
ics  (Givoli  and  Keller,  1990),  (Givoli,  1992).  These  boundaries  pertain  to  a  two-dimensional  circular 
geometry,  and  would  be  difficult  to  implement  in  a  general  manner.  Kallivokas  (Kallivokas  et  ai,  1991) 
have  developed  a  high-frequency  (singly  asymptotic)  non-reflecting  boundary  for  two-dimensional  ap¬ 
plications  that  is  based  on  elastodynamic  ray  theory  (Keller,  1958). 

Doubly  asymptotic  approximations  were  originally  developed  to  treat  structures  submerged  in  acoustic 
media  by  Geers  in  the  1970’s  (Geers,  1971),  (Geers,  1978).  Since  that  time,  DAA’s  for  acoustic  me¬ 
dia  has  been  extended  and  more  systematically  derived  (Felippa,  1980a),  (Geers  and  Felippa,  1983), 
(Nicholeis-Vullierme,  1991),  (Geers  and  Zhang,  1991).  A  first-order  DAA  has  also  been  formulated 
for  electromagnetic  scattering  by  Geers  and  Zhang  (Geers  and  Zhang,  1988).  The  first-order  doubly 
asymptotic  approximation  for  infinite  elastic  media  was  formulated  heuristically  by  Underwood  and 
Geers  (Underwood  and  Geers,  1981),  and  further  developed  by  Mathews  and 
Geers  (Mathews  and  Geers,  1987).  As  mentioned  previously,  doubly  asymptotic  approximations  are 
much  more  robust  than  singly  asymptotic  approximations. 
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SECTION  3 


FUNDAMENTAL  EQUATIONS 


In  this  section  we  present  the  fundamental  equations  that  are  used  to  formulate  the  DAA’s.  Small  char¬ 
acters  with  arrow  accents  represent  vectors,  capital  characters  represent  tensors  or  matrices,  and  stan¬ 
dard  font  characters  represent  scalars. 

Let  us  express  the  displacement  vector  u{x,t)  through  a  Helmholtz  decomposition  in  terms  of  a  scalar 
potential  <i>{x,t)  and  a  vector  potential  V»(x,y,z,t)  as  (see,  e.g.,  (Achenbach,  1973)), 

u{x^t)  =  V(j){x,t) +  V  X V-xf{x,t)  —  0.  (3.1) 


The  wave  equation  for  a  uniform  elastic  medium  then  separates  into  the  uncoupled  wave  equations 

CdVV(x,<)  = 


(3.2) 


4vV(^,0  =  tpix.t). 


in  which  cq  and  cs  are  the  dilatational-  and  shear-wave  speeds,  respectively,  given  by 


2  "I"  2^  2  ^ 

C£)  —  ,  Cg  —  — , 

P  P 

where  A  and  fi  are  the  Lame’  parameters  for  the  medium.  The  stress  tensor  is  given  by 
(see,  e.g.,  (Eringen  and  Suhubi,  1975)) 


(3.3) 


a  =  -1-  2// VV(^  -I-  VV  x  ^ -|-  V  x  ■ifV). 


(3.4) 
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Where  I  is  the  identity  tensor. 


An  exact  integral-equation  solution  to  (3.1)  and  (3.2)  on  the  smooth  surface  or  surfaces  of  an  elas¬ 
tic  medium  is  provided  by  Love’s  integral  identity,  which  may  be  written  in  Laplace-transform  space 
as  (Cruse  and  Rizzo,  1968) 

y{F,s)  +  J^ii{P,s)f{F,P,s)dSp  =  Jj{P,s)U{F,P,s)dSp.  (3.5) 


where  the  usual  Cartesian  tensor  notation  applies,  R  is  defined  as  |  fp  —  r/r  |,  and  where  ^(i?,  s)  and 
xiR,  s)  are  given  for  three  dimensions  as: 
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Figure  3- 1 .  Description  of  semi-infinite  domain  for  force  normal  to  surface. 


The  second-order  tensor  operators  used  in  equation  (3.5)  for  a  half-space  have  been  derived  previously 
by  Banerjee  and  Mamoon  (Banerjee  and  Mamoon,  1990).  The  first  group  of  operators  is  based  on  a 
force  vector  normal  to  the  half-space  boundary,  and  the  second  group  is  based  on  a  force  vector  parallel 
to  the  boundary.  Figure  3-1  shows  the  coordinate  system  and  the  other  geometric  quantities  used  for 
the  first  group  of  operators  ,  which  obviously  pertain  to  an  axisymtric  problem.  The  notation  used  is 
that  of  Banerjee  and  Mamoon. 
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Z  S  0  Plane 


Figure  3-2.  Description  of  semi-infinite  domain  for  force  parallel  to  surface. 

The  operator  components  pertain  to  the  following  point  forces: 

1.  Single  dynamic  force  at  (0,0,+c) 

2.  Single  dynamic  force  at  (0,0,-c) 

3.  Dynamic  double  forces  at  (0,0,-c) 

4.  Dynamic  centre  of  compression  at  (0,0,-c) 

5.  Line  of  dynamic  centers  of  compression  extending  from  2  =  — ct0  2  =  — oo 

6.  Dynamic  doublet  at  (0,0,-c) 

The  resulting  operator  components  are  given  the  appendix. 

Figure  3-2  shows  the  coordinate  system  for  the  second  group  of  operators.  The  point  forces  consist  of 


1 .  Single  dynamic  force  at  (0,0,+c) 

2.  Single  dynamic  force  at  (0,0,-c) 
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3.  Dynamic  double  force  with  moment  at  {0,0,-c) 

4.  Line  of  dynamic  double  forces  with  moment  extending  from  z  =  -c  to  z  -  -oc 

5.  Dynamic  doublets  at  (0,0,-c) 

6.  Line  of  dynamic  doublets  extending  from  z  =  —c  to  z  =  — oo  with  strength  proportional  to  the 
distance  from  z  =  —c 

The  resulting  operator  components  are  given  in  the  appendix. 

The  components  in  (A.l)  and  (A.2)  are  identical  to  those  of  the  first  of  (3.6),  which  is  the  standard 
Kelvin  solution  for  a  point  force  in  an  infinite  elastic  space.  The  other  terms  are  additional  terms  to  be 
added  to  the  basic  Kelvin  solution  to  form  the  half-space  solution  (Brebbia  and  others,  1984).  The  T 
operators  for  the  half-space  are  found  from  (A.l)  and  (A.2)  by  differentiating  the  U  components  with 
respect  to  spatial  coordinates  (Brebbia  and  others,  1984).  The  operators  (A.l),  based  on  the  cylindrical 
coordinate  system,  are  transformed  to  the  Cartesian  coordinate  system  using  standard  tensor  techniques 
for  coordinate  transformations. 
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SECTION  4 


EARLY-TIME  APPROXIMATIONS 


Early-time  approximations  are  non-reflecting  boundary  relations  valid  for  ct  a,  where  a  is  a  char¬ 
acteristic  dimension  of  the  boundary  S  and  c  is  a  wave  speed  (Felippa,  1980b).  The  approximations 
are  local',  they  depend  only  on  geometrical  information  in  the  neighborhood  of  the  field  point  F.  Thus, 
early-time  approximations  for  an  infinite  domain  (whole-space)  are  identical  to  those  for  a  semi-infinite 
domain  (half-space). 

We  present  herein  two  derivations  of  early-time  approximations.  The  first  derivation  is  based  on  ray 
elastodynamics  (Keller  (Keller,  1958)).  The  second  is  based  on  work  by  Felippa  (Felippa,  1980b)foran 
acoustic  medium,  which  is  in  turn  based  on  Kirchhoff’s  retarded  potential  formula  (Kirchhoff,  1883). 
The  two  derivations  produce  the  same  relationships  involving  the  displacement  potentials.  Those  re¬ 
lationships  are  then  used  to  formulate  two  early-time  approximations  (ETA’s)  involving  tractions  and 
displacements. 

4. 1  RAY  ELASTODYNAMICS  DERIVATION. 

In  this  section  we  present  the  first  method  for  deriving  the  early-time  approximations.  This  derivation  is 
accomplished  using  ray  elastodynamic  theory  developed  by  J.B  Keller  and  H.B.  Keller  (Keller,  1964), 
(Ahluwalia  et  al,  1969),  (Achenbach  and  others,  1982).  The  basic  assumption  of  ray  elastodynamics 
is  that  disturbances  are  propagated  along  rays  to  which  the  disturbance  surface,  i.e.  the  wavefront,  is 
always  normal.  Figure  4-1  shows  the  basic  geometry. 
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Ray  Bundle 


Figure  4- 1 .  Definition  of  ray  bundle  for  elastodynamics. 


The  potentials  in  (3. 1)  may  be  represented  by  a  ray  series  (Ahluwalia  et  ai,  1969) 


<^(f,5)  =  exp[-sv?(f)/cD]0^^o(  — ) 


=  exp[-5v?(f)cs]0^^o(— ) 

cs 


(4.1) 


where  the  phase  function  (f{x)  and  the  amplitude  functions  ^>’"(f)  and  ^'"(f)  depend  on  the  spatial 
variables,  but  are  independent  of  s,  the  Laplace  transform  parameter.  Keller  (Keller,  1964)  has  shown 
that,  at  high  frequencies,  the  displacement  potentials  can  be  adequately  described  with  only  the  first 
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term  in  each  expansion  (4. 1 )  as 


(p{F,Ti,s)  «  exp[-.s///cD] 
xI){F,t},s)  th  exp[-S7;/c5] 


[(7/  +  R^){r)  +  R(;) 

R(^  R(^ 


,  i 

2 


[(7/  +  R^){r]  +  Rc)\ 


(i>{F,0,s) 

i’{F,0,s), 


(4.2) 


where  the  argument  0  denotes  the  initial  wavefront  on  the  surface,  R^  and  are  the  principal  radii 
of  curvature  of  the  initial  surface,  and  r]  is  the  coordinate  along  the  ray,  which  is  orthogonal  to  each 
wavefront  and  straight  for  homogenous  media.  Equation  4.2  relates  the  displacement  potentials  at  some 
point  T)  along  the  ray  to  the  initial  disturbance  surface,  and  to  the  initial  amplitude  functions  ^(F,  0,  s) 
and 


At  a  field  point  F  on  the  surface,  the  ray  is  directed  along  the  normal  to  the  surface.  Therefore,  taking 
derivatives  of  0  with  respect  to  the  local  coordinates  at  F,  we  obtain 

R^R(^ 


^{F,r],s)  ^  — ^exp[-577/cD] 

On  CD 


n  i 

2 


1 


-  -exp[-sr]/cD] 


[(77  +  R^){tj  +  R(i) 

R^R(;(2rj  +  +  R^) 

{tj  +  R^Y{7]  +  R^y 


(j){F,0,s) 


(7/  +  Ri)iri  +  Rc) 

R^Rq 


n  i 


(^(F,0,5), 


d(l> 

if  ==  '’■ 

#  «  0. 
Si 


(4.3) 


From  the  first  of  these,  the  partial  derivative  of  <j)  with  respect  to  the  local  normal  at  the  point  F  on  the 
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surface  is  just 


^(F,0,s)  =  lirn^(F,7/,5)  ^  )</>(  F,  0,  s ) 

an  n-^0  an  cq  IR^Kc^ 


—  ~( - 1- «)i?^(F,  0,  s), 

CD 

where  k  is  the  mean  curvature  (positive  when  the  surface  is  convex).  In  (4.4),  the  derivative  of  the  scalar 
potential  at  point  F  is  related  to  the  scalar  potential  itself  at  the  same  point  F.  We  can  readily  obtain  a 
similar  result  for  the  vector  potential  0.  Thus  we  have  the  following  results  for  the  partial  derivatives 
of  the  displacement  potentials  with  respect  to  the  local  coordinates  at  point  F: 


-{ —  +  K)<i>{F,s) 
CD 


-{-  +  K)rl^{F,s) 
Cs 


:{F,s) 


These  results  are  used  below  to  form  early- time  approximations. 
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tangent  plane  at  F 


Figure  4-2.  Coordinate  systems  on  surface  for  early-time  approximations. 

4.2  RETARDED  POTENTIAL  DERIVATION. 

In  this  section  we  develop  the  elastodynamic  equivalents  to  Felippa’s  first  two  early-time  approxima¬ 
tions  for  acoustic  media  (Felippa,  1980b).  We  use  two  local  coordinate  systems,  as  seen  in  Figure  4-2. 
The  first  is  a  Cartesian  system  centered  at  point  F,  where  x  and  y  define  the  plane  tangent  to  the  sur¬ 
face  at  F.  At  point  F,  x  and  y  coincide  with  the  orthogonal  surface  coordinates  and  (.  The  second 
coordinate  system  is  a  cylindrical  system  (r,  6,  c)  with  =  a:^  -t-  y^. 

The  distance  between  the  field  point  F  and  the  integration  point  P  is  given  by  R,  which  isR  = 

Felippa  expanded  the  c  coordinate  of  the  surface  in  a  Taylor  series  about  the  point  F  to  get 

c  =  AR^  +  BR^  +  {C-  A^)R-*  +  {D-  ^A^B)R!>  4-  •  •  •  .  (4.6) 
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where 


A  =  7—^  cos^  0  +  —■  sin  0  cos  9  +  -r-r  sin^  0 
Ox^  uxuy  oy^ 


^3  -w 

B  =  7-77  cos^  0  +  377 — 7 — 7-  sin  6  cos^  9  +  ■  ■  ■  -—77  sin^  6 
ax'^  axaxuy  Oy^ 


The  derivatives  . . .  are  assumed  to  exist  at  F,  i.e.,  the  surface  is  smooth. 

First,  we  discuss  irrotational  motion.  The  integral  solution  to  the  wave  equation  for  the  scalar  displace¬ 
ment  potential  4>  is 


27r<?i(F,0  =  J{R 


_^d(b{PAR)  ,  ry-2^R 


+  dSp, 


which  is  Kirchhoff’s  retarded  potential  formula  (Love,  1904a),  (Eringen  and  Suhubi,  1975),  wheref/?  = 
t  —  represents  the  retarded  time. 


We  expand  0  about  the  field  point  F  in  a  Taylor  series  to  get 


«PI)  =  4,(F,t)  +  ^(F,t)i+^{F,t)(+l-^{F,t}e 


I  (p  .X  cr  .  L^(f  t)  r2  4. 


At  F,  the  coordinates  ^  and  (  coincide  with  the  x  and  y  coordinates,  and  can  be  approximated  as 


e  =  7?cos0[l  +  0(H'*)],  C  =  Rsme[l  +  0{R^)] 


(4.10) 
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Therefore, 


(/.(P,0  =  cf>{F,t)  +  ^^{F,t)  Rcose[l  +  0(R')]  +  ^i^{F,t)  Rsme[l  +  0{R')] 

ax  ay 

+  I \^^{F,t)  cos^  B  +  ^(FJ)  cos0sin0  (4.11) 

yiOx^  ox  ay 

The  normal  derivative  of  (j)  may  be  similarly  expanded  to  obtain 

+  ^^{F.t)Rsin0[l  +  O{R^)] 
anoy 

We  now  introduce  (4.11)  and  (4.12)  into  (4.8)  and  utilize  Felippa’s  expansions  for  R~  US  and  R-^^dS, 
viz. 


R-US  =  {1  +  ]^AlR^  +  ---)dRde 


(4.13) 


R-^^^dS 

an 


-{A  +  2BR  +  3(C  -  A^)R'^  +  ■■■  )dRd6 
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to  obtain 


rCDt  r'iir  f)0 

27r(^(F,0  sa  /  /  —{FjR)dOdR 

Jo  Jo  on 


(4.14) 


+ r‘  r  {<f>iF,tR)+—^{F,tR)]{-A)dedR. 

Jo  Jo  I  CD  ) 

Note  that  (?i>(F,  in)  and  its  normal  derivative  are  functions  only  of  the  retarded  time  and  not  the  co¬ 
ordinates  R  and  0. 


Next,  we  integrate  over  0  and  divide  through  by  to  get 

<p(F,0  «  j^'''^{FM)dR-Kj^'''\^cf>iF,tR)  +  ^kF,tR)]dR  (4.15) 

where  k  is  again  the  mean  curvature  of  the  surface  at  F.  Then,  making  the  variable  change,  R  =  CD{t- 
tR),^^  get 


(f>{F,t)  ^  CD  ■^{F,iR)dtR  —  cdk  i^(l>{F,tR)  +  {t  —  Ir) 


(4.16) 


Finally,  integrating  by  parts  and  noting  that  (f>{F,  0)  =  0  for  quiescent  initial  conditions,  we  obtain 


^{F,t)  -I-  KCD(l>{F,t)  «  -CD-^{F,t).  (4.17) 

which,  when  transformed  to  the  Laplace  transform  domain,  gives 

{s  +  kcd)(I){F,s)  ^ -cd^{F,s).  (4.18) 


We  observe  that  this  is  identical  to  the  first  of  (4.5) 


19 


Now  we  consider  solenoidal  (rotational)  motion.  The  integral  solution  to  the  wave  equation  for  the 


vector  potential  0  is 


(4.19) 


+R-^ 


m 

dn 


R  ' 

+  -HP.tR) 
cs  y 


dS. 


This  is  again  Kirchhoff’s  retarded  potential  formula  (Eringen  and  Suhubi,  1975).  Thus,  application  to 
i/’  of  the  procedure  described  between  (4.8)  and  (4.18)  yields 


(5  +  KCs)tf{F,s)  =  -cs-^{F,s). 


(4.20) 


We  observe  that  this  is  identical  to  the  second  of  (4.5). 


4.3  FORMULATION  OF  ETA2 . 

In  this  section  we  employ  (4.5)  to  obtain  early-time  relations  between  surface  tractions  and  displace¬ 
ments  that  do  not  involve  spatial  derivatives.  For  (4.5)  we  used  a  local  curvilinear  coordinate  system 
C,  n  with  the  origin  at  point  F,  characterized  by  two  principal  curvatures,  R^  and  Rc^  (Figure  4-1).  Thus 
we  can  write  ^  =  R^^Oi  and(  =  72^02,  so  that  our  metric  coefficients  are  A 1  =  72^, /i2  =  72^,and/i3  =  1. 
We  also  note  that,  because  the  spatial  derivatives  of  the  scalar  and  vector  potentials  are  non-negligible 
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only  for  the  normal  direction,  we  have  the  early-time  surface  formulas  (Moon  and  Spencer,  1961) 


V(?i>  « 

on 


V^j/’  =  V(V  ■  Ip)  —  'V  X  V  X  tjj 


+lS  +  2-=^  +  iRc  -  RdRl^R^’i’S 


a^4,  84, 

dn^  dn 


(4.21) 


0  0  0 


VV4> 


0  0  0 


V(VxV')+(Vxi/>)Vs 


<»  <»  ^ 
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We  obtain  the  early-time  displacement-potential  surface  relations  by  introducing  the  first  two  of  (4.21) 
into  (3.1),  which  yields 


(4.22) 


We  obtain  the  early-time  traction-potential  surface  relations  by  introducing  the  last  three 

of  (4.21)  into  (3.4)  and  defining  surface  tractions  and  as  the  negatives  of  cr„^,  and  cr„^  on 

the  surface;  this  yields 


i 

Cn 


,  2  «  .d<i>. 


(4.23) 


Also,  we  obtain  early-time  surface  relations  between  the  potentials  and  their  normal  derivatives  by  in¬ 
troducing  the  third  and  fourth  of  (4.21)  into  (3.2)  and  applying  the  Laplace  transform  for  quiescent 
initial  conditions  to  get 


^  “'"an 


(4.24) 


[^  +  2.?^  +  {R,-R0RfRrU 


With  (4.5)  and  (4.22)  -  (4.24)  in  place,  we  may  now  proceed  to  derive  the  desired  early-time  relations. 
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First,  we  eliminate  0  in  the  first  (4.5)  by  introducing  the  first  of  (4.24),  eliminate  in  the  resulting 
equation  by  introducing  the  first  of  (4.23),  eliminate  in  the  resulting  equation  by  introducing  the  first 
of  (4.22),  and  make  use  of  the  first  of  (3.3);  this  yields 


(S  +  HCD)tn  ~  {pCDS  +  4//Ac)5U„, 


(4.25) 


where  terms  of  0(5°)  on  the  right  side  have  been  dropped. 

We  find  the  early-time  stress-displacement  relations  involving  the  shear  stresses  by  first  using  the  sec¬ 
ond  of  (4.5)  to  eliminate  ^  and  ^  in  the  last  two  of  (4.24),  which  yields 


dn? 


4-  2k.{scs^  +  k)  ~  ^(7?^  ^  -  R\)]rp^ 


(4.26) 


dn^ 


«  [.s^C5^  +  2k(scs^  +  k)  —  R^^{R^  ^ 


Second,  we  use  the  second  of  (4.5)  again  to  eliminate  ^  and  ^  in  the  last  two  of  (4.23),  which  yields 


D-l 


U  ~  -^c  )V’c] 


(4.27) 


02  / 


Third,  we  eliminate  and  in  (4.27)  by  introducing  (4.26)  to  get 


23 


^  p4 '  +  {2k  -  R^^)scs^  +2{ K^  -R^\ R^ '  -  R^')] 

-R'^^{k  +  R^  ‘)]i/’(: 

(4.28) 

«  -pclis^cs^  +  {2k  -  R-^^)scs^  +  2{k^  -  R^\R^^  -  R^')} 

-R-\K  +  R-^)\4’i. 


Fourth,  we  use  the  second  of  (4.5)  a  third  time  to  eliminate  ^  and  ^  in  the  last  two  of  (4.22),  which 
yields 


RS  (SC5*  +  K —  /?(. ')?/’< 


(4.29) 


«  —{scs^  +  K  —  R^^)ip^. 


Fifth,  we  eliminate  and  in  (4.28)  by  introducing  (4.29)  to  obtain 

{s  +  l3cs)t^  «  pcs{s^  +  Rj^css  +  l3R^^)u^ 


(4.30) 


{s-^cs)t^  ^  pcs{s‘^  +  R(-^css  -  l3R^^)u^, 
where  j3  =  k  —  R^^  =  —{k  —  Rj^). 

The  traction-displacement  relations  (4.30)  are  not  quite  in  proper  form.  Thus,  we  multiply  the  first 
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of  (4.30)  by  +  /ics)  and  the  second  by  -  fics)  to  get 

(3  +  K.cs)t^  ^  pcs ^  ^  (3^  +  R7^Css)U(^ 

3  +  pcs 


(4.31) 


{s  +  KCs)t^  ^  pcs^ +  Rr  ^Css)u^, 
s  -  pcs 

where  terms  of  0(3°)  on  the  right  side  have  been  dropped.  Finally,  we  employ  long  division  and  retain 
terms  only  in  3^  and  3  to  obtain 


(3  +  KCs)^^  ~  [/)Cs3  +  2k//]3U^ 


(4.32) 


(s  +  KCs)t(i  W  [pcs5  +  2kp]sU(i 


The  desired  relations  (4.25)  and  (4.32)  may  be  assembled  to  produce  the  second-order  early-time  ap¬ 
proximation  ETA2  for  an  elastic  medium; 

(3  -(-  KCo)tn  =  {pCDS  +  4/iK:)3U„ 


(s  +  KCs)t^  =  {pcss  +  2pK)su^ 


(4.33) 


(3  -I-  KCs)t^  =  {pcss  -f  2pk)sU(^ 


These  are  immediately  inverse-transformed  to  obtain  the  ETA2  in  the  time  domain: 


tn  4-  KCotn  —  pCoUn  +  ApKUn 


4-  KCst^  =  pcsu^  4-  2pKu^ 


(4.34) 


4-  KCst(^  =  pcsu^  4-  2pKU^ 
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which  may  be  written 


t{F,t)  +  KCt{F,t)  —  pCu(F,t)  +  KAu{F,t) 


(4.35) 


where 


C  = 


CD 


cs 


cs 


,  A  = 


ip 


2p 


(4.36) 


For  very  early  times  (5  — ^  oo)  the  second  order  ETA  obviously  reduces  to  the  first-order  ETA 


t{F,t)  =  pCu{F,t) 


(4.37) 


We  have  completed  the  formulation  of  the  early-time  (high-frequency)  approximations  ETAi  and  ETA2. 
These  approximations  are  commonly  called  the  plane  wave  and  curved  wave  approximations,  respec¬ 
tively.  A  key  property  of  these  approximations  is  that  the  coefficient  matrices  are  diagonal;  therefore, 
the  approximations  are  spatially  local  with  uncoupled  local  components. 

We  now  rotate  the  early-time  approximations  from  local  coordinates  to  global  coordinates.  Displace¬ 
ments  and  tractions  in  local  coordinates  can  be  rotated  into  global  coordinates  as 

u{F,t)  =  Q{F)u{F,t) 


t{F,t)  =  QiF)i{F,t), 

where  Q  is  the  rotation  tensor.  Then  (4.37)  becomes 

Qt{F,t)  =  pCQu[F,t), 


(4.38) 


(4.39) 
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which  yields  ETAi  in  global  coordinates: 

i{F,t)  =  pQ  (4.40) 

where,  because  Q"*  =  Q^,Cl=  Q^CQ.  The  corresponding  first-order  high-frequency  approximation 
is  therefore 

t{F,s)=pastliF,s),  (4.41) 

The  second-order  early-time  approximation  can  be  similarly  rotated  to  get  ETA2  in  global  coordinates: 

LiF,t)  +  KQi{F,t)  =  pQuiF.t)  +  K\ii{F,t),  (4.42) 

where  A  =  Q^^Q-  Hence  the  second-order  high-frequency  approximation  is 

{s  +  /cC)i(F,  5)  =  {pQs  -h  kK)su{F,  $)  (4.43) 

From  this  point  we  will  work  in  global  coordinates  only;  the  underlines  will  be  dropped  with  the  un¬ 
derstanding  that  all  quantities  are  global. 
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SECTION  5 


LATE-TIME  APPROXIMATIONS  FOR  ELASTIC  DOMAINS 


Late-time  approximations  pertain  to  low-frequency  response,  for  which  the  chacteristic  elastodynamic 
wavelengths  c©//  and  cs//  are  much  larger  than  the  characteristic  dimension  of  the  boundary.  We  will 
derive  first  the  late  time  approximations  (LTA’s)  for  the  infinite  domain  (whole-space)  and  then  derive 
the  LTA’s  for  the  semi-infinite  domain  (half-space). 


5.1  LATE-TIME  APPROXIMATIONS  FOR  A  WHOLE-SPACE. 


We  use  (3.5)  -  (3.7)  to  construct  an  approximation  that  approaches  the  exact  solution  as  5  — >  0.  First, 
we  expand  (3.7)  in  Taylor  series  as 


7(7?,  s)  = 


1  Cc  I  Co  2  7?  c|  3  2  3\ 

(tt  +  1)  “  0(73"  +  )’ 


27? 


D 


rd)  cs 


8'c' 


D 


Cs 


(5.1) 


R,  d 


1 


D 


8 


Note  that  the  second  of  these  has  no  0(5^)  term,  and  that  the  0(5^)  term  in  the  first  is  independent  of 
R.  Substitution  of  these  equations  into  (3.6)  produces 


U{F,P,s)  =  U^{F,P)  +  sU\F,P)  +  Ois^)  (5.2) 

f{F,P,s)  =  f^{F,P)  +  0{s^). 
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Where  the  components  of  U^,  U^,  and  are  given  by 


U'  -  1  [1-2,/  ^(l-Wl 

''  I27r,i(l-u)  [  2cd  cs  J  " 

=  ~3:r(l  -  v)B?  ~  ~  ~ 

With  these  relations,  we  can  now  form  the  late-time  approximations  for  the  whole-space. 


5.1.1  First-Order  Late-Time  Approximation:  LTA^ . 


Let  us  introduce  (5.2)  into  (3.5)  and  retain  terms  only  of  0(s°).  Transformation  back  into  the  time 
domain  then  yields  the  quasi-static  relation 

iu(F,0  +  j^u{P,t)f\F,P)dSp  =  Jj{P,t)U^{F,P)dSp,  (5.4) 

which  we  recognize  as  the  standard  Somigliana  identity  for  elastostatics.  The  relation  (5.4)  is  based  on 
the  Green’s  function  for  an  infinite  elastic  medium  obtained  by  Kelvin  in  1848  (Thomson,  1948).  With 
the  spatial-operator  definitions 


Mpt)  =  f  q{P,t)U\F,P)dSp 

Js 


tq[P,t)  =  f  q{P,tmF-F)  +  f<‘{F,P)]dSp, 

J  s 

where  6(F  —  F)  is  the  Dirac  delta-function  and  '  denotes  a  spatial  tensor-integral  operator,  (5.4)  may 
be  expressed  as  LTA^^ : 


t{P,t)  =  F-^fu(F,f). 
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LTA|'^  is  a  spatialy  non-local,  singly  asymptotic,  quasi-static  approximation. 


5.1.2  Second-Order  Late-Time  Approximation:  LTA^- 

Now  we  introduce  (5.2)  into  (3.5)  and  retain  terms  through  O(s^).  This  yields 


i{?(F,  s)  +  u{P,  s)f°(F,  P)dSp  =  ^  t{P,s)U^iF,  P)clSp  +  5  ^  t{P,  s)U^iF,  P)dSp. 

(5.7) 

With  the  spatial-operator  definitions  (5.5)  and  the  new  spatial-operator 

Aq{P,t)=  l^qiP,t)U\F,P)dSp.  (5.8) 

(5.7)  may  be  expressed  as  the  low-frequency  approximation  LFA“: 

{sA  +  B)t{P,  s)  =  f  u(F,  5).  (5.9) 

Transformation  of  this  equation  back  into  the  time  domain  produces  the  second-order  Late-Time  Ap¬ 
proximation  for  the  whole  space  LTA^ :. 


At{P,  t)  4-  Bt{P,  t)  =  fu(F,  t) 


(5.10) 


LTA^  is  a  spatialy  non-local,  singly  asymptotic  approximation. 

5.2  LATE-TIME  APPROXIMATIONS  FOR  A  HALF-SPACE. 


In  this  section  we  derive  the  late-time  approximations  for  a  semi-infinite  domain,  or  half-space.  We 
will  use  the  same  method  as  that  used  in  the  previous  section  for  the  whole-space. 

Only  the  A,  B,C,  P,  Q,  and  S  coefficients  given  in  the  appendix  are  dependent  on  s.  We  need  to  expand 
these  coefficients  in  Taylor  series  about  5  =  0  to  find  the  first-  and  second-order  late-time  approxima¬ 
tions.  The  expansions  for  the  unique  terms  are: 
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Ai  — 


cs^  \  cs^r: 

'2R\  '  2cd^R\)  '  \3cs  2cd^J 


1  ^  C. 


2  \  ^_2 
+ 


+  0(s") 


As  — 


Bs  = 


Cs  = 

B,  = 

^6  = 
Be  = 
Ss  = 


-3  ' 


Cs 


0R2^  '  ■2cd^R2' 


-  s^  +  0{s^) 


bCs^  8  CD 


^  ^  2R2^^  2cd^  R2V  '  I8C52  SCD 

-1  3c5^  \  ,  /  -1  3c5^ 


2R2^  ^  2cd^  R2V  '  VScs^  8  CD 


+ 


-  s^  +  0{s^) 


-1 


cs 


2R2^  2cd^R2^)  V8cs2  Scd^ 

-1  3cs2 


3  3cs2 

+ 


+ 


2R2^  2  cd^R2^J  \Scs^R2  8cd^R2. 


+  0(s") 


-1 


Cs 


0R2^  2cd^R2\ 

'  -I  2  ' 

_!_  +  _« _ 

.2X2^  2cd^R2\ 


+ 


+ 


'  +  +  0(^=) 


8  cs^  i?2  8  cp'*  R2 


-3  cs' 


8cs^  8cd‘‘ 


5'  +  0(/) 


(5.11) 


The  terms  of  0(s°)  are  those  that  produce  the  first-order  LTA  for  an  elastic  half-space.  The  expres¬ 
sion  for  Ai  contains  the  only  O(s^)  term.  Thus,  only  U^,,  Ul^,  Ul^,  and  contribute  to  the  second- 
order  LTA.  We  note  that  the  tensor  operator  f{F,P)  involves  spatial  deriviatives  of  the  components 
of  U  (Brebbia  and  others,  1984),  (Cruse  and  Rizzo,  1968).  Additionally,  the  expression  for  Ai  is  in¬ 
dependent  of  position  (no  R  dependency).  Therefore,  f{F,  P)  will  contribute  only  0(s°)  terms  to  the 
second-order  LTA,  as  is  the  case  for  the  whole-space  formulation. 


31 


5.2.1  First-Order  Late-Time  Approximation:  LTA(^. 

The  first-order  LTA  for  a  half-space  is  a  quasi-static  approximation  expressed  by  (5.4),  but  with  T°{F,P) 
and  U°{F,  P)  given  by  expressions  corresponding  to  the  half-space  Green’s  function  of 
Mindlin  (Mindlin,  1936),  rather  than  that  of  Kelvin  (Brebbia  and  others,  1984).  These  expressions  are 
constructed  by  augmenting  (5.3),  the  Kelvin  solution,  with  additional  terms  (Brebbia  and  others,  1984). 
The  additional  terms,  for  both  the  U°{F,P)  and  T°(F,  P),  are  the  0(s°)  terms  of  IP  through  (7®, 
[see  (A.l),  (A.2)]. 

For  U^{F,P),  the  augmentation  is 

UfiiF,  P)  =  f7°(F,  P)  4-  F*(F,  P)  (5.12) 
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in  which  the  components  of  U*{F,  P)  are; 


where  A'<f 


u:.  = 


U’  = 


u:,  = 


u-  =  K, 


8(1  -//)^-(3-4;y)  + 


(3  —  Au)r\  —  2cx  6ca;rj 


m 


+ 


.,A 


(3  —  4i/)ri  4(1  —  i/)(l  —  '2i/)  ^  6ca:ri 


^  ^  (3-4i/)r|  2cx  ( 1  -  -irl ' 


Rl 


m  V  ^2 


4r(l  —  u){\  —  2u) 


1  - 


=  Rdr^rz 

U‘ 

rzu:. 

-9 

r2 

M  ^ 

CO 

II 

r2 

u:. 

=  U’ 

yx 

^y 

Rly 

=  Ad  1 

r  +  ri  y  R2{r  +  ri) 

(3  —  4/2)  4(1  — /2)(1  —  2//)  6cx 


Rl 


Riir  +  ri)2 


(3  —  4//)r3  2cx 


Rl 


+ 


Rl 


1  - 


Rl 


4r(l  —  //){!  —  2i>) 
r  +  ri 


1  - 


^2(r  +  ri)/J 


(5.13) 


iMl^andCsee  Figure  5-1) 


n  =  z{F)  -  z{P) 
r[  =  z{F)  -  z{P') 
R  =  (r,Ti)5 
X  =  z{F) 


r2  =  x{F)  -  x{P) 
r',  =  xiF)  -  x{P') 
R2  =  (r'r')5 
c  =  z{P) 


^3  =  y{F)  -  y{P) 
r'  =  y{F)  -  u{P') 
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Figure  5- 1 .  Three-dimensional  geometry  (in  case  of  the  half-space,  the  infinite  free  surface  lies  in  the 
X  -  y  plane). 

T{F,P)  is  augmented  as 

THiFP)  =  f{F,P)  +  f-{F,P)  (5.14) 

in  which  the  components  of  T*{F,P)  are  given  by 

T‘^  =  a’kiTik  2  =  1,  X  =  2,  y  =  3  (5.15) 

where  rik  is  the  surface  normal  (defined  as  positive  going  into  the  medium),  and  where  the  components 
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of  (7*  are,  with  Kt  = 


(7.,,  = 


O’,,,,  = 


zyz 


yyz 


(7, 


yx 


Ill 

Ri 

Rtvi 

Ri 

Ktr3 
m  L 


1  —  2i/ri  — 


l-2u- 


1  -  2i/  - 


3(3  —  4i/)aT'|  —  3crj(5.r  —  c)  SOcxr'^ 


Rl 


Ri 


2  '*^2 
3(3  —  4i/).Trj  —  9cx  —  3c^  30cxr'j 


Ri 


Ri 


3(3  —  4i/).tr',  —  9car  —  3c^  30ca;r’'^ 


Rl 


Ri 


Ill 

Rl 


(1  —  2t')(3ri  —  \ur\)  — 


3(3  —  4j2)rir2 

m 


+6crj 


,  (1  —  ‘2ux  —  2i/c)  30cxr|rj 


Rl 

4r^(l  —  i/)(l  —  2u) 


Ri 


r  +  ri 


1  - 


i?2(r  +  n)  RIJ\ 


xyz 


Ktr2rz 

Rl 


3(3-4z/)ri  A{\~u){\-2v)  (  1 


Rl 


+ 


1  \  30cxr( 
+  -7r]-  ' 


Rl  L 


(1  —  2j/)(3ri  —  4t/rj) 


r  +  ri  Vr  +  Ti  '  i?2/  -^2 

3(3  —  AuYirl  6crj[(l  —  2i')x  —  2vc\ 


Rl 


Rl 


“idcxT^r^  4r^(l  —  j/)(l  —  2i/) 


Ri 


r  +  ri 


1  - 


R2{r  +  rj)  RIJ\ 


A>2 

m  L 


3(3  —  4i/)r?  6c  /  \  /  bxr'l 

-{l-2u)--^^-L±  +  —  U+{l-2u)r[+  ^ 


Rt  L  o  ^  3(3-4j/)r|7;  6c  f_,  ,  2 

=  (l-2i/)ri - ^2 - -  (1  -  - 


Ri  . 

Sxr^Tj 

Rl 


I<tr2r3 

R\ 

A>2 


—3(3  —  Au)r[  +  6c  (  1  —  2z/  + 


bxr[ 

~W 


Rl 


3- 


(1  -2i/)(5-4j/) 


r^(3A2  +  r'l)' 
Ri{r  +  ri)  , 


3(3  —  4i/)r2  4r’^(l  —  i/)(l  —  2i/) 


Rl 


(r  +  ri)2 


(3,  -  3  -  2.ri  +  tri) 


(5.16) 
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V(r3  r  3(3-4t/)r^  Ar^{\  -  i'){\  - '2u)  (  rj(3i?2  +  r\)\ 

Hi  [  m  {r  +  ny  V  mir  +  n)  ) 


6cx  I  orj 

■  m  I  Ri 


{I -2u){3-Au)  - 


3(3  -  4//)r2 


4r^(l  -  ;/)(l  -2u)  /  _  r|(3^2_+^\ 
(r  +  ri)2  \  Ri{r  +  ri)  ) 

+^(c-(l-2.K  +  |r?) 


^zzy  -  ^112 

^2 


^zxy  —  <^132 


<yy  = 


Z{Z-Au)rlr\  _  ^  .  2  _ 

Rl  Rl  ^  ^  Rl 


(1  -2i/)(3-4i/)  - 


3(3  -  4//)r2 


4r^(l-^/)(l-2z/)  /  _  rl{3R2  +  r[)\ 
(r  +  ri)2  Ri(r  +  rj)  j 


+  -S  fc-(l  -2//)r;  + 


A>2  [  3(3  -  4iy)ri  -  j/)(1  -  2u)  f  rl{3R2  +  r'jX 

Ri  [  Ri  (r  +  ri)2  \  Ri{r  +  r,)  ) 


6cx  /  5ro 


(1  -2//)(5-4i/)  - 


3(3-4j/)r|  4r2(l  -  z/)(l  -  2i/)  /  r^(3A2  +  ri) 


(r  +  ri)^ 


Ri{r  +  ri) 


+  ^  -  (3  -  2z/)r'i  +  . 


In  operator  notation,  LTAi  for  the  semi-infinite  domain  appears  the  same  as  that  for  the  infinite  domain, 
given  by  (5.6). 
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5.2.2  Second-Order  Late-Time  Approximation:  LTA.2^. 

The  second-order  late-time  approximation  includes  both  0(5°)  and  0(5^ )  terms  in  the  half-space  op¬ 
erators.  In  (A.l)  and  (A.2),  the  components  pertain  to  the  Kelvin  solution.  The  U'^  components, 
‘  which  are  the  image  complements  to  contribute  components  that  have  0(5^)  terms.  From  the  dis¬ 

cussion  following  (5.11),  there  are  no  0(s* )  components  of  the  T  operator.  Thus,  the  operators  for  the 
second-order  approximation  are  [c.f.,  (5.2)] 

Uh{F,P)  =  Uff{F,P)  +  Ulj{F,P)  +  Ois^) 


(5.17) 


f„iF,P)  =  fff{F,P)  +  Oisf. 


When  we  substitute  the  results  of  (5.17)  into  (3.5),  we  get 


~it{F.s)  +  l^u{P.s)f°{F,P]dSp  =  Jj{P,s)0^„{F,P)dSp  +  s  lj{P,s)U},iF,P)dSp 


(5.18) 


We  can  express  this  equation  using  the  spatial-operators  definitions  of  (5.5)  and  (5.8)  as  the  second- 
order  low-frequency  approximation  for  a  half-space 

LFA^  :  [sA  +  B]t{P,s)  =  fu(F,5).  (5.19) 

Inverse-transforming  to  the  time  domain  we  obtain  the  second-order  late-time  approximation  for  the 
half-space. 


LTA^  :  At{P,t)  +  m(P.  t)  =  ru{F,s).  (5.20) 

We  have  now  completed  the  formulation  of  both  the  first-order  and  second-order  late-time  approxima¬ 
tions  for  the  half-space.  A  notable  result  is  that  the  late-time  approximations  for  the  half-space  have 
the  same  forms  as  their  counterparts  for  the  whole-space. 
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SECTION  6 


DOUBLY  ASYMPTOTIC  APPROXIMATIONS 


In  this  chapter  we  combine,  in  systematic  fashion,  the  early-time  and  late-time  approximations  into 
first-order  and  second-order  doubly  asymptotic  approximations  for  both  whole-  and  half-spaces.  We 
conclude  with  the  discretization  of  the  DAA  relationships  using  the  boundary  element  method. 

6. 1  FIRST-ORDER  DOUBLY  ASYMPTOTIC  APPROXIMATIONS. 

Because  the  only  differences  between  the  first-order  early-  and  late-time  approximations  for  infinite 
and  semi-infinite  domains  reside  in  the  operators  f{F,P)  and  6^  ( F,  P ) ,  we  can  formally  develop  first- 
order  DAA’s  for  the  two  domains  simultaneously.  We  use  the  method  of  operator  matching  for  this 
purpose  (Nicholas- Vullierme,  1991),  (Geers  and  Zhang,  1991).  The  appropriate  trial  equation  is 

t{P,s)  =  [sih  +  Uo]u[F,s)  (6.1) 

where  the  unkown  spatial  operators  Uq  and  Ux  are  not  functions  of  s. 

For  s— 0,  we  write  (6.1)  as 

t{P,s)  =  \Uo  +  0{s)]u{F,s)  (6.2) 

and  match  with  (5.6)  as  s— >  0,  which  yields 

(Jo  =  B-^t.  (6.3) 

This  is  the  asymptotic  match  for  the  static  limit. 

For  s-^  cx),  we  write  (6.1)  as 


t{P,s)  =  [III  +  0{s  ^)]su{F,s) 

(6.4) 

and  match  with  (4.41)  as  s— 

^  oo  to  get 

Ux  =  pC. 

(6.5) 

Introducing  this  result  and  (6.3)  into  (6.1),  we  obtain,  in  transform  space. 

DAAxis)  : 

f(P,5)  =  [pCs  +  B-'^t 

u{F,s) 

(6.6) 

and  in  the  time  domain 

DAAxit)  : 

t{P,t)  =  pCu{F,t)  -b  B~ 

^tu{F,t) 

(6.7) 

This  result  was  heuristically  formulated  by  Underwood  and  Geers  (Underwood  and  Geers,  1981).  Note 
that  the  DAAi  for  elastic  domains  is  not  spatially  local,  because  of  the  second  term  on  the  right,  which 
is  the  result  of  5  — >  0  matching. 
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6.2  SECOND-ORDER  DOUBLY  ASYMPTOTIC  APPROXIMATIONS. 


In  this  section  we  develop  second-order  DAA’s  consisting  of  a  two-term  low-frequency  component  and 
a  two-term  high-frequency  component.  For  the  reason  given  in  the  previous  section,  we  can  develop 
the  whole-space  and  half-space  DAA’s  simultaneously..  The  appropriate  trial  equation  is 

{s  +  fo)t{P,s)  =  (s^pC  -h  sU,  +  Uo)u{F,s) 

in  which  highest-frequency  matching  has  already  been  performed,  and  To,  f/i,  and  Uo  are  unknown 
spatial  operators  that  are  not  functions  of  s. 

We  first  match  the  trial  equation  as  s  — >  0  with  the  low  frequency  approximations  (5.9)  and  (5.19).  The 
left  side  of  these  may  be  inverted  to  yield 

t{P,s)  =  [I  -sB-'^ A +  0{s^)]B-^tuiF,s),  (6.8) 

inasmuch  as  {B  4-  =  (^  +  sB~^  A)~^ B~^  =  [^  ~  sB~^A  +  0{s^)]B~^.  Similarly,  the  left  side 

of  (6.8)  may  be  inverted  to  obtain 

t{P.  s)  =  [I-  sfo^  +  0{s'^)]f^AUo  +  SU,  +  3^pC)u{F,  s).  (6.9) 

We  can  first  match  (6.8)  and  (6.9)  to  0(s°)  to  obtain 

fo^Uo  =  B-^t.  (6.10) 


and  then  match  to  0(s^ )  to  get 


-  To-'Uo)  =  -b-^Ab-^t. 


(6.11) 


We  now  match  the  trial  equation  as  s  — >■  cxj  to  the  second-order  high-frequency  approximation,  (4.43), 
which  we  write  as 

[I -H  5“^/c(7]f(F,  s)  =  {pC  +  s~^kA.)su{F,s)  (6.12) 

Note  that  we  have  already  matched  at  0(s°).  To  match  at  0(s“^),  we  divide  (6.8)  through  by  to  get 

{1  +  s~^to)t{P,s)  =  {pC  +  s~^Ui  +  s~'^Uo)su{P,s)  (6.13) 

The  left  sides  of  (6.12)  and  (6.13)  may  be  inverted  as  5  ^  oo  to  get 

t(F.s)  =  [I  —  s~^kC  +  0{s~^)][pC  +  s~^kA]su{F,s) 

t{F,s)  =  [I-s-^fo  +  0{s-^)][pC  +  s-'^Ui  +  0{s-^)]su{F,s)  (6.14) 
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Thus,  matching  the  0(s  ‘ )  terms,  we  obtain 


{/i  _  pfoC  =  kA  -  pkC^ 


(6.15) 


From  (6.10),  (6.1 1)  and  (6.15),  the  three  unknown  spatial  operators  may  be  found  as 


To  =  [5-*f  +  KipC^  -  A)]  {pC  +  =  n 


Ui  =  B-^t  -  ClB-^ AB-^t  (6.16) 

Uo  =  nB-^t. 

Substitution  of  these  results  back  into  the  trial  equation  produces  the  second-order  DAA  in  s-space: 


DAA2{s):  [s  +  n]t  =  [s^pC  +  s{B-^t-hB-^AB-^t)  +  nB-^t]u.  (6.17) 


Inverse  transformation  then  yields  the  second-order  DAA  in  the  time  domain: 

DAA^it):  't  +  nt  =  pC{i+{B-^t-ClB-'^AB-^t)il+ClB-^tu.  (6.18) 


6.3  MIXED-ORDER  DAA. 

Here,  a  one-term  low-frequency  matching  is  combined  with  a  two-term  high-frequency  matching  to 
form  the  mixed  doubly  asymptotic  approximation  DAAi_2. 

The  trial  equation  is  again  (6.8).  For  the  low-frequency  match  we  use  (5.9)  with  A  =  0.  The  matching 
procedure  of  the  previous  section  then  yields  (6.16)  with  A  =  0.  Thus  we  obtain  the  mixed-order  DAA 


DAA,.2{s)  ; 

[s -I-  ^1-2]^  =  -f  S.S  ^f-t-J7i_2.S  ^f]u, 

(6.19) 

DAAi_2(0  : 

t  +  Cli-2t  =  pCu{t)  +  B~^tu{t)  +  Q\-2B~^tu{t). 

(6.20) 

where  17i_2  = 

p-^B-^tC-^  +  k{C  -  ^-^AC-I)]. 

6.4  BOUNDARY  ELEMENT  DISCRETIZATION. 

A  robust  method  for  applying  DAA’s  to  general  surfaces  is  finite  element  discretization 
methods  (Underwood  and  Geers,  1981),  (Mathews  and  Geers,  1987).  We  will  discretize  the  ETA2  first, 
by  invoking  the  Bubnov-Galerkin 
method  (Cook  et  aL,  1989).  We  have,  then 

u{F,t)  =  N(F)u(f) 

(6.21) 

t{F,t)  =  N{F)t{t). 
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where  a  boldface  quantity  denotes  a  computational  (n-dimensional)  vector  or  matrix.  In  (6.21),  u(i)  is 
the  vector  for  nodal  displacements  of  the  entire  discretized  surface,  t(<)  is  the  corresponding  vector  of 
nodal  tractions,  and  N(F)  is  a  matrix  of  shape  functions  that  interpolate  nodal  values  from  the  nodes 
to  the  point  F.  Applying  Bubnov-Galerkin  discretization  to  (4.42),  we  obtain  the  matrix  ETA2 

t(0  +  Ft(0  =  Du(0  +  Lu(0,  (6.22) 

where 

F  =  J  N^(F)kCN(F)  dS 
D  =  J  N^(F)pCN(F)  dS 
L  =  J  N^(F)/cAN(F)  dS, 

Similarly,  the  matrix  ETAi  follows  from  (4.40)  as 

t{t)  =  Du(f).  (6.23) 

Discretization  of  the  late-time  approximations  proceeds  in  the  same  way,  so  that  (5.10)  yields  the  matrix 

LTA2 

Ai{t)  -f  Bt(0  =  Gu{t)  (6.24) 

where  A,  B,  and  G  are  the  discretized  forms  of  (5.8)  and  (5.5).  The  matix  LTAi  is  obtained  by  taking 
A  =  0  in  this  equation. 

Matrix  doubly  asymptotic  approximations  may  be  obtained  from  (6.22),  (6.23),  and  (6.24)  by  the  method 
of  matrix  matching,  which  proceeds  in  exactly  the  manner  as  that  described  in  Chapter  6  for  operator 
matching.  The  resulting  matrix  DAA2  is 


t(0  Wt(0  =  Dii(0  +  (B-*G  -  WB-*  AB-'G)u(0  -f  WB-^Gu(0, 


(6.25) 


where 


W  = 


B-‘G-FD-L 


D  +  B-'AB-^G 


-1 


The  DAAi_2  is  found  from  (6.25)  by  setting  the  matrix  operator  A  to  zero.  The  DAAi  is  found  by 
matching  only  the  0(s°)  terms,  it  is 


t(0  =  Du(0  +  B-^Gu(0  (6.26) 

Note  that  we  can  premultiply  the  nodal  traction  vector  t  ( < )  by  the  diagonal  surface-area  matrix  contain¬ 
ing  element  areas  to  get  the  nodal  force  vector. 

In  this  section  we  derived  the  DAAi,  DAA2,  and  mixed-order  DAAi_2  matrix  formulations  used  on 
the  operator  matching  methods  presented  in  section  6.2.  The  matrix  formulations  pertain  to  both  the 
whole-space  and  half-space.  The  first-order  matrix  DAA  is  a  first-order  ordinary  differential  equation; 
the  second-order  DAA  and  the  mixed-order  DAA  are  second-order  ordinary  differential  equations. 
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SECTION  7 


MODAL  SOLUTIONS  FOR  AN  ELASTIC  WHOLE-SPACE 


In  this  section  we  develop  exact  n  =  0  and  n  —  1  modal  temporal  impedance  relations  for  a  spherical 
cavity  in  an  elastic  whole-space  and  the  corresponding  DAA2  relations.  We  show  that  the  second-order 
modal  DAA  agrees  with  the  exact  temporal  impedance  relation  for  n  =  0. 


7. 1  EXACT  MODAL  EQUATIONS  FOR  n  =0. 


This  first  modal  derivation  is  based  on  the  acoustic  work  of  Geers  and  Zhang  1991 
(Geers  and  Zhang,  1991).  We  start  with  a  spherical  surface  in  an  elastic  whole-space  and  assume  mo¬ 
tion  in  the  radial  direction  only.  The  following  equations  apply: 


Uo  = 


d(i>o 

dr 


cro 


(A +  2^) 


dr'^ 


-f 


2A  d(l>o 
r  dr 


(7.1) 


<i>Q  _  d^<f>Q  2  d(l>o 
Cq  dr"^  r  dr 

Here,  (po  is  the  scalar  displacement  potential  for  the  n  =  0  mode,  uo  and  oq  represent  the  radial  dis¬ 
placement  and  normal  traction,  respectively. 


Laplace  transforming  the  third  of  (7.1)  and  multiplying  through  by  r^,  we  obtain 


2^ Vo  5<^o 

r - h  2r - 

dr^  dr 


-(-)Vo  =  0. 

CD 


(7.2) 


This  equation  is  a  form  of  Bessel’s  equation  with  77  =  0  (Arfken,  1970);  the  solution  is 


Mc,s)  =  Ms)ko(—)  ^  fo{s)ko{0,  (7.3) 

CD 

where  ^  /o(.s)  is  an  unknown  function  of  the  Laplace-transform  variable,  and  A:  ((f)  is  the  mod¬ 

ified  spherical  Bessel  function  of  the  zeroth  order. 

Using  the  third  of  (7.1)  to  eliminate  from  the  second  of  (7.1),  Laplace  transforming,  and  then  in¬ 
troducing  (7.3),  we  get 


ao{r,s)  =  ps'^fo(s)ko{0 - -Uo 


(7.4) 


Now 


dkoU) 

dr 


(7.5) 
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(7.6) 


where  thus,  (7.3)  and  the  first  of  (7.1)  yield 


tio(r,s)  =  Ms)  — 

CD 


From  (7.4)and  (7.6),  we  obtain 


Mc.s)  =  -  —]uo{r,s), 

where  the  prime  denotes  a  derivative  with  respect  to  But  (Arfken,  1970), 

koio  = 

Kio  =  -c-Hr'  +  r^) 


koiO 

k'oU) 


i  +  r 


so  that  (7.7)  may  be  expressed  as 


/  ,  CD  ,  2  ,  4/15  4/iCD  , 

(s  +  — )po  =  ipCDS  + - [-  — 5— )uo,  (7.9) 

where  the  traction  po  =  —  o’o  and  4  has  been  replaced  by  This  is  the  exact  modal  impedance 
relation  forn  =  0  motion  of  a  spherical  surface  in  an  elastic  whole-space. 


Note  that,  at  early  times  (s  —>■  00),  (7.9)  yields 


Po  =  PCDUo, 


(7.10) 


which  is  the  standard  plane- wave  temporal  impedance  relation  (Love,  1904a).  At  late  times  (s  -+  0), 
(7.9)  yields 


4p 

Po  =  — , 
r 


(7.11) 


which  is  the  static  solution. 


We  will  now  develop  an  n  =  0  modal  impedance  relation  for  shear  traction  and  corresponding  tangen¬ 
tial  displacement.  Starting  with  a  spherical  surface  in  an  elastic  whole-space  and  assuming  uniform 
tangential  motion  that  is  only  a  function  of  the  radial  spatial  coordinate,  the  vector  displacement  be¬ 
comes  a  scalar  function.  Therefore,  the  following  equations  apply; 


Vq  = 


_  d'^tpo  2  dipo 
dr^  r  dr  ’ 

^  2dxpo 

c%  dr"^  r  dr 


(7.12) 
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where  7/’o  is  the  scalar  displacement  potential  for  the  tangential  n  =  0  motion,  and  vo  and  tq  represent  the 
tangential  displacement  and  the  shear  stress  on  the  surface,  respectively,  for  the  n  =  0  mode.  Laplace 
transforming  the  third  equation  of  (7.12)  and  multiplying  through  by  r^,  we  obtain 


r  -TT  +  - (  — )  i'o  =  0- 

ar^  or  cs 


(7.13) 


This  equation  is  a  form  of  Bessel’s  equation  with  n  —  0.  Following  the  same  process  as  that  used  for 
radial  n  =  0  motion,  we  obtain  for  tangential  n  =  0  motion 


Cs. 


2fis 


{s  H - )7'o  —  [/>C5S  + 

V  V 


H - ^h’o- 


(7.14) 


We  have  now  derived  the  n  =  0  modal-impedance  relations  for  both  radial  and  tangential  motions.  In 
the  next  section  we  show  that  DAA,  when  adapted  to  these  one-dimensional  problems,  gives  identical 
equations. 


7.2  DAA2  EQUATIONS  FOR  n=0. 

We  now  develop  a  modal  second-order  doubly  asymptotic  approximation  (DAA^ )  for  an  elastic  whole- 
space.  Again,  we  follow  the  acoustic  work  of  Geers  &  Zhang  (Geers  and  Zhang,  1991),  outlining  the 
process  for  the  radial  motion  case  only,  as  the  tangential  case  follows  identical  steps. 

A  second-order  high-frequency  approximation  (HFA2)  is  obtained  by  the  introduction  of  the  last  of  (7.8), 
into  (7.7),  which  yields 


HFA2  :  {s  +  ^)po  =  {pCDs'^  +  (7.15) 

This  equation  is  identical  to  the  exact  n  =  0  modal  equation  because  all  the  ray-theory  assumptions  are 
satisfied  exactly  for  n  =  0  motion. 

The  second-order  low-frequency  approximation  (LFA2)  is  found  from  the  expansion  of  the  last  of  (7.8) 
as  5  — >  0,  which  produces 

=  =  (7.16) 

When  we  substitute  this  result  into  (7.7)  with  ,  and  keep  terms  of  (9(s^ )  and  0(s°),  we  get 

LFA2  :  <7'o(r,5)  =  -(O5  -I-  —)iio(r,s).  (7.17) 

r 

Because  there  are  no  terms  of  0(.s^),  the  second-order  low-frequency  approximation  (LFA2)  for  the 
n  =  0  modal  solution  is  identical  to  the  first-order  low-frequency  approximation. 

To  form  the  modal  second-order  DAA  for  an  elastic  whole-space,  a  trial  equation  is  selected  as 

(5  -f  7o)po  =  {pCD^'^  +  U\S  -f-  Uo)uo,  (7.18) 
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where  Tq,  U\,  and  Uq  are  unknown  constants  that  are  not  functions  of  s,  and  po  is  the  pressure,  i.e., 
(To  =  — Po.  Solving  for  po  in  equation  (7.18)  we  get 


Po  —  (.s  +  To)  +  l-hs  +  Uo)uo. 


(7.19) 


Now  as  .s  0, 


(s  +  To)  ^  —  To  ^  —  To  +  To  +  •  •  •  . 


(7.20) 


Substituting  this  result  into  (7.19),  keeping  terms  up  to  and  including  0(5^),  and  matching  the  terms 
with  (7.17),  we  get 


T„-'Uo  = 

r 


Ui-To^Uo  =  0 


To  match  the  high  frequency  terms,  we  arrange  the  trial  equation  (7.18)  for  5  — *•  oo  as 


(7.21) 


Po  =  (1  -  To5  ^  +  To  5  ^  )[pCD  +  UiS  ^  +  Uo$  ^]5Uo. 


(7.22) 


A  similar  rearrangement  of  (7.15)  yields 


_  Cqs  ^  /  C£)  ^2  —2  \r  I  —1  I  4pC£)  _2.| 

Po  —  (1 - 1-  ( — )  s  ...  )[pc£)  H - s  q - - — s  ]suo. 


(7.23) 


Matching  (7.22)  and  (7.23)  through  order  3  ^  we  obtain 


Ui  -  pcDio  = - p — 

r  r 


(7.24) 


Simultaneous  solution  of  (7.21)  and  (7.24)  yields 


rr 

-to  —  — 

r 


Uo  = 


4pc|) 


(7.25) 


u,  = 

r 


so  that  (7.18)  gives  for  radial  n  =  0  motion. 


DAA2~° -radial 


(3  -1-  — )po  =  (pCD^^  +  —S  +  ^^^)Uo, 


(7.26) 
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Note  that  this  equation  is  identical  to  the  exact  equation  (7.9).  Proceeding  in  the  same  manner  for  tan¬ 
gential  n  =  0  motion,  we  find 

DAAJ^^-tangential  (s  -| — -)to  =  {pcss^  -f  — .s  -f  — )t’o,  {J .'ll) 

r  r 

which  is  identical  to  (7.14). 

Regarding  the  assembly  (7.26)  and  (7.27)  as  one  doubly  asymptotic  approximation  for  n  =0  motion 
of  a  spherical  surface,  we  conclude  that  DAA^"”  is  an  exact  relation. 


7.3  EXACT  MODAL  EQUATIONS  FOR  n=  1 . 


Now  we  will  derive  the  n  =  1  modal  solution  of  a  spherical  cavity  in  an  infinite  elastic  medium  us¬ 
ing  direct  transformation.  We  start  with  relationships  between  stresses  and  potentials  for  axisymmetric 
motions  (Eringen  and  Suhubi,  1975). 


cr  =  -|-  2p 


d^(j)  ^ 

dr'i  Qr 


T  = 


‘  dr  do 


2^ 

r  do 


d 

do 


d^{rip) 

dr^ 

d'^irip) 

dr'^ 


a.2S) 


1  d'^{rtp)  d 

r  dOdr  dr 


1  di^irxl}) 
r  dOdr 


The  relationships  between  the  radial  and  meridional  displacements  and  the  displacement  potentials  are: 


u 


V 


dr 


d‘^{rxl^) 

dr^ 


r'V'^'tp 


I  d(i>  1  d^irtp) 
r  do  r  dOdr 


(7.29) 


We  can  represent  the  displacement  potentials  for  n=  1  as 

<t>  =  /l(5)^•x(OPl(cos0)  (7.30) 

^  =  5'i(5)/:i(77)Fi(cos^), 

where  ^  ~  A:iisa  modified  spherical  Bessel  function.  Pi  is  Legendre  polynominal  and 

/i  and  gi  are  coefficients  dependent  only  on  the  Laplace-transform  variable  s.  Multiplying  the  above 
equations  by  Pm{cosO),  integrating  over  the  surface,  and  invoking  the  orthogonality  property  of  Leg¬ 
endre  functions,  we  get,  after  some  rearrangement  and  substitution; 


(Ti  —  AVj<?!)i  -f  2/i 


d‘^4>i 

dr^ 


dr 


d^iri’i)  „2 


dr^ 


-  rVitpn 


(7.31) 


n  =  - 


r  I  or  r 


d^iri^i)  __„2 


dr"^ 


—  rV^ipi 


,  ld{ripi)  d 

- 'p - 

r  dr  dr 


1  d{ripi) 


dr 
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where 


*  dr'^  r  dr  r'^ 


(7.32) 


We  can  now  eliminate  the  second-order  partials  by  using  the  wave  equations  and  Vj?/’!  = 


2ji  4// 

(Ti  =  /)s  - III - t’l  (7.33) 

r  r 

2  ,  2/^  2/z 

ri  =  ps^i’i - ui - til, 

r  r 

We  can  get  the  n  =  1  relationships  for  the  displacements  (7.30)  by  employing  Legendre  functions 
and  (7.32). 


til  =  -x - 'H’l 

or  r 


dil^i  t/)i  <1)1 

Vi  —  -z — I - , 

or  r  r 


We  can  use  (7.29)  and  the  above  equations  to  solve  for  /i  and  gi  as 


/i  = 


9\  = 


-cpri-cskijri)  -  rsk[{i]))ui  -f  2rcDCski{r})vi 

-2cDCski{T])ki{^)  +  rscski{r])k'i{()  -|-  r^s^k'i{rj)k'i{^) 

_ rcocski{^)ui  +  r'^scsk[{^)vi _ 

-2cDCski(T])ki{^)  -I-  rscski{r])k'i{^)  +  r^s^k'i{r})k'i{^)' 


(7.34) 


(7.35) 


We  can  eliminate  ©i  and  t/^i  in  (7.33)  using  /i  and  and  (7.30).  Substitution  of  the  expansions  for  the 
modified  spherical  Bessel  functions,  given  by 

h(o  =  e-'(r'+r‘) 
i'i(«)  =  e-'(-2r"-2r"-r‘) 
hM  =  +’)'') 


yields  the  traction-displacement  relations. 


,2  Ti  To,  f  CTi  \  ^  3  ,  U2  2  ,  Ui ,  Uo,  f  til 


+  :^S  +  ^] 


—  [/tC^  -| - s  -| — —s  -(- 


j.  J.3 


(7.36) 
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where 


Ti 

To 

U2 

Ui 

Uo 


2cd  +  cs  0 
0  2cd  +  Cs 


=  (2co  +  C5)I 


2c^ 


D  + 
0 


0 

2c^  +  c 


(24  +  4)1 


4//  +  p(?o  +  /’Cdcs  4/i  -  'Ipcocs  1 

2p  -  pc  DCs  2p  +  /)c|  +  2pcDCs  J 

Ap[2cD  +  C5]  +  pieces  +  cdc|]  Ap[2cD  +  cs]  -  2p[cpCs  +  cdc|] 
2p[2cD  +  Cs]  -  p[cdCs  +  cd4]  ‘^p[^cd  +  cs]  +  2p[cl)Cs  +  cd4] 

4424  +  4]  +  444]  4/z[24  +  4]  -  2444]  1 
2424  +  41-444]  2/.[24  +  41  +  2444]  I 


(7.37) 


Equation  7.36  is  the  exact  n  =  1  modal  equation.  We  will  now  derive  various  approximate  relations 
for  n  =  1  motion. 


7.3.1  n=  1  High-Frequency  Approximations. 

First,  we  write  (7.36)  as 

[I  +  r-'Ti5-'  +  r-2To5-2]ti  =  [pC  +  r-'Ujs"'  +  +  r-^Uo^'^l^Uj, 

(7.38) 

where  ti  and  Ui  represent  the  n  =  1  traction  and  displacement  vectors.  To  get  HFA2“\  we  drop  terms 
of  0(s“^),  and  0(s~^)  inside  both  brackets;  to  get  HFAi=\  we  drop  terms  of  0(s“^),  0(s“^),  and 
0(s“^)  inside  both  brackets.  Thus  we  have 

HFAr'  :  [I^  +  -Ti]ti  =  [pCs^  +  — ]ui,  (7.39) 

r  r 

HFA^=i  :  ti  =  pCsUi  (7.40) 

7.3.2  n=l  Low-Frequency  Approximations. 

We  premultiply  (7.38)  through  by  ^  =  r^l{2c\,  +  C5)  to  get 

[1  +  -^5  +  (  — Iti  =  (rco)~4Uo  +  ?’Ui5  +  r^U2S^  +  r^pCs^jui  (7.41) 

Cl  Co 

where  Cq  =  2c£,  +  c\  and  ci  =  2c£,  +  c|/2c£»  +  cs.  Now  we  divide  through  by  1  +  to 

get,  through  O(s^), 

LFA^=^  :  ti  =  [Ki  +  K'lsjui  (7.42) 
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(7.43) 


where,  with  =  (cdcs/cq)^  and  cl  =  cxCdcs{cd  +  cs)lcl-, 

4//  +  pc\ 

2/i  -  pci  2p  +  2pcl 


K,  =  {rc^)-uo=-^ 


if  4//  +  pci  4p  -  2pcl 


K',  =  cr'lco'^CiU,  -  cj'Uol  =  /.cr'(c^  -  4)  I  ./  I  (7.44) 

If  we  keep  terms  only  through  0(s°),  we  get  the  first-order  low-frequency  approximation 

LFA^^  :  ti=KiUi,  (7.45) 

7.4  n=  1  DOUBLY  ASYMPTOTIC  APPROXIMATIONS. 

As  the  form  of  the  modal  first-order  approximations  for  low-  and  high-frequencies  are  similar  to  the 
form  of  the  first-order  low-  and  high-  frequency  general  DAA  (section  6.1),  we  can  use  the  same  process 
to  combine  (7.40)  and  (7.45)  to  get  the  first-order  DAA  for  the  n  =  1  motion. 

DAAx~^  =  /oCui  -f  KiUi  (7.46) 

We  can  form  a  one-low,  two-high  DAA  (DAA”f])  from  (7.39),  (7.45)  and  the  trial  equation 

(5  -t-  To)ti  =  (5^pC  +  sUi  -f  Uo)ui 

Proceeding  in  the  manner  described  in  section  6.2,  we  get  for  the  spatial  operators 

To  =  [^^-25  +  K,11pC1->  =  W, 

r  r 

Ui  =  Ki  (7.47) 

Uo  =  WiKi, 

which  yields 

DAA^zl  {s  -f  Wi)ti  =  {s^pC  -H  sKx  -h  WiKi)ui  (7.48) 

For  the  full  second-order  DAA  (DAA5=' ),  we  again  start  with  the  trial  equation  (7.47)  and  employ  (7 .39) 
and  (7.42).  Using  the  matching  procedures  employed  to  obtain  (7.21),  we  find  from  low-frequency 
matching 

To^Uo  =  Ki  (7.49) 

To'Ui-To"Uo  =  K'l. 
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The  first  equation  is  from  the  0(s°)  match,  and  the  second  equation  is  from  the  0(5^ )  match. 

For  high-frequency  matching,  we  use  the  procedure  employed  to  obtain  (7.24)  this  yields 

U, -TopC  =  r-‘U2-  r-'TipC  (7.50) 

which  results  from  the  (s~^ )  match.  Finally,  simultaneous  solution  of  (7.49)  and  (7.50)  yields 

To  =  [5^  -  —  +  K,1[,>C  -  K',]-' =  V, 

r  r 

Ui  =  (7.51) 

r  r 

Uo  =  ViKi. 

which  yields  DAA2  for  n  =  1  motion 

DAAr'  (5  +  Vi)tx=  6 -f  ViKi  uj  (7.52) 

In  this  section  we  have  developed  exact  and  DAA  temporal  impedance  relations  for  n  =  0  and  n  —  1 
motions  of  a  spherical  surface  embedded  in  an  infinite  elastic  medium. 
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SECTION  8 


CANONICAL  PROBLEMS 


Here  we  compare  DAA-based  and  exact  results  for  canonical  problems,  as  done  previously  by  Geers 
and  Underwood  (Underwood  and  Geers,  198 1)  and  by  Mathews  and  Geers  (Mathews  and  Geers,  1987). 
In  a  recent  analytical  study,  Oberai  (Oberai,  1994)  generated  DAA-based  response  for  step-wave-excited 
infinite-cylindrical  and  spherical  shells  embedded  in  infinite  elastic  media.  His  comparisons  with  cor¬ 
responding  exact  results  showed  DAAi  to  be  fairly  satisfactory  and  DAA2  to  highly  satisfactory. 

We  first  discuss  the  computer  programs  developed  to  implement  the  DAA  equations  and  then  we  com¬ 
pare  results  for  four  problems: 

1 .  A  spherical  cavity  embedded  in  an  elastic  whole-space  subjected  to  an  n  =  0  internal  step  pres¬ 
sure. 

2.  A  spherical  cavity  embedded  in  an  elastic  half-space  subjected  to  an  n  =  0  internal  step  pressure. 

3.  A  spherical  cavity  embedded  in  an  elastic  whole-space  subjected  to  an  n  =  1  internal  step  trac¬ 
tion. 

4.  A  uniform  step  pressure  applied  to  a  portion  of  the  surface  of  a  half-space. 

8.1  DAA  COMPUTER  PROGRAMS. 

We  developed  three  computer  programs  for  computing  DAA-based  responses.  The  first  program, 
EDAA 1 ,  implements  the  theory  of  first-order  doubly  asymptotic  approximations  for  three-dimensional 
elastic  media,  both  infinite  and  semi-infinite.  The  second  program,  EDAA  1-2,  implements  the  theory 
of  the  mixed-order  (one-low,  two  high)  doubly  asymptotic  approximation  for  three-dimensional  infinite 
elastic  media.  The  third  program,  EDAA2,  implements  the  theory  of  second-order  doubly  asymptotic 
approximations  for  three-dimensional  elastic  media,  again  for  both  infinite  and  semi-infinite  spaces. 
The  routines  for  the  first-order  kernels,  U°  and  T”  (6.24),  are  based  on  software  obtained  from 
Mathews  (Mathews,  1990). 

Eight-node  serendipity  elements  are  used  to  discretize  the  surface  (Hughes,  1987)(  Figure  8-1).  The 
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Figure  8-1.  Eight  node  serendipity  element. 


weighting  functions  are 

^^2(6, 6) 

iV4(6,6) 
NsUui2) 
NeUu^2) 
-^7(61 6) 
iV8(6,6) 

where  ^1  G  [-1, 1]  and ^2  €  [-1, 1]. 


—  —^(1  “  6)(1  “  6)(1  +  6  +  6) 

=  j(l  +  6)(l  —  6)(“i  +  6  +  6) 

=  ^(l+6)(l-e2)  (8.1) 

=  ^(1  +  6)(i  +  6)(“i  +  6  +  6) 

=  ^(i-e?)(i  +  6) 

=  ^(1  —  6)(i  +  6)(— 1  -  ^1  +  6) 

=  ^(i-6)(i-e2) 

The  methods  explained  in  Section  6.4  are  used  to  form  the  kernel 


matrices  with  3  x  3  Gauss  quadrature  employed  for  element  integration.  The  required  matrix  inversions 
are  performed  with  standard  Gauss  elimination  techniques.  Due  to  the  basic  assumption  of  linearity, 
all  matrices  are  formed  only  once,  and  decomposed  once,  to  form  the  DAA  equations.  The  loads  f 
and  f  are  calculated  at  each  step.  The  programs  use  trapezoidal  time  integration  techniques.  The  first- 
order  DAA  (DAAi)  is  a  first-order  ordinary  differential  equation  for  u  and  the  mixed-order  (DAAi_2) 
and  second-order  DAA  (DAA2)  second-order  ordinary  differential  equation  for  u.  Therefore  the  two 
programs  have  different  implementations  of  the  trapezoidal  method. 

The  DAAi  (EDAAl)  matrix  equations  (6.26)  are  advanced  in  time  to  step  n  +  \  using  the  following 
algorithm  (Hughes,  1987); 

St . 

U„+1  =  u„  -h  — u„ 


D  +  — B  G  1  Un+l  — 


fn+1  + 


^n+l  —  ^(Un+1  —  U„+i), 

where  St  is  the  time  interval  used,  and  Un+i  is  the  predictor. 

The  programs  EDAAl -2  and  EDAA2  use  the  trapezoidal  method  for  second-order  differential  equa¬ 
tions.  The  trapezoidal  method  that  is  used  for  time  integration  in  EDAA2  is  shown  below: 


4d  +  |:(B-'G  -  WB-*AB-'G)-^WB-ig1u„+i 

.0t‘‘  ot 


=  (Wf  -f-  f)  +  (B-'G  -  WB-' AB-*G)(^u„  +  u„) 

vt 

4  4 

-f-  D(— u„  +  — u„  4- u„) 
ot^  ot 


^n+l  —  ^n) 


U71+1  —  Un 


The  time  integration  method  for  EDAAl -2  is  formed  from  the  EDAA2  by  dropping  terms  containing 
A. 

8.2  SPHERICAL  CAVITYAVHOLE-S PACE/STEP  PRESSURE  PROBLEM. 


The  first  problem  is  a  spherical  cavity  embedded  in  an  infinite  elastic  medium  and  excited  by  an  in¬ 
ternal  step  pressure.  This  problem  possesses  radial  symmetry,  and  has  a  well-known  analytical  solu¬ 
tion  (Timoshenko  and  Goodier,  1970).  With  a  as  the  cavity  radius  and  po  as  the  pressure  magnitude. 
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the  radial  displacement  of  the  cavity  is  given  by 


where 


u[t)  =  — “‘(cos +  —  sin 

4/i  (  /i 

— ae~“‘(—  sin  afit  +  cos  a3t) 


1 


+  1  —  e  cos  a dt  +  —  sin  a 3t) 


Cpjl  -  2i/) 
a(l  —  i/) 


The  corresponding  analytical  DAAi  solution  is  simply 

UDAAit)  =  ^  {l  -  . 


(8.4) 


(8.5) 


(8.6) 


The  boundary-element  model  of  the  cavity  boundary  consisted  of  24  eight-node  elements  over  the  en¬ 
tire  spherical  surface.  The  analytical  exact,  analytical  DAAi ,  and  numerical  DAAi  solutions  are  shown 
in  Figure  8-2  for  the  parameters  p  —  1.00,/i  =  1/6,  r/  =  1/4,  a  =  1,  and  po  =  1  •  The  analytical  DAAi 
and  numerical  DAAi  solutions  eire  seen  to  be  almost  identical,  and  the  DAAi  solutions  agree  well  with 
the  analytical  exact  solution  at  both  early  and  late  times.  As  previously  observed  by  Underwood  eind 
Geers  (Underwood  and  Geers,  1981),  the  DAAi  solutions  do  not  exhibit  the  response  overshoot  seen 
in  the  exact  solution.  This  is  expected  from  first-order  differential  equations  like  (6.7). 

The  radial  displacement  response  from  the  second-order  DAA2  (EDAA2)  is  compared  with  the  ex¬ 
act  response  in  Figure  8-3.  Also  in  this  figure  is  the  DAAi  boundary-element  result.  The  agreement 
between  the  second-order  DAA  and  the  exact  analytical  solution  is  excellent.  The  small  difference  be¬ 
tween  the  two  results  is  due  to  element  discretization,  which  yields  a  discrepancy  in  surface  area  of  1%. 
The  second-order  DAA2  yields  the  “overshoot”  characteristic  of  elastic  media  problems. 

Figure  8-4  shows  a  comparison  of  exact,  boundary  element  DAAi,  mixed-order  modal  DAAi_2,  and 
mixed-order  boundary  element  DAA  (EDAA12)  results,  for  the  cavity  subjected  to  an  internal  step 
pressure.  The  exact  Timoshenko  and  the  modal  DAAi_2  results  are  identical.  This  is  because  LTAi  is 
identical  to  LTA2  for  this  case. 


8.3  SPHERICAL  CAVITY/  WHOLE-SPACE/n=  1  STEP  TRACTIONS  PROBLEM. 

This  problem  involves  a  spherical  cavity  embedded  in  an  infinite  elastic  medium  and  excited  by  an 
n  =  1  internal  step  traction.  The  problem  is  axisymmetric.  The  exact  modal  relation,  (7.36),  includes 
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Time  CQt/a 

Figure  8-2.  Exact  and  DAAi  results  for  spherical  cavity  in  an  elastic  medium  subjected  to  an  internal 
step  pressure. 


Figure  8-3.  Exact,  DAAi,  and  DAA2  radial  displacement  of  a  spherical  cavity  in  an  infinite  elastic 
medium  subjected  to  an  internal  step  pressure. 
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Figure  8-4.  Exact,  DAAi,  and  DAAi_2  radial  displacement  of  a  spherical  cavity  in  an  infinite  elastic 
medium  subjected  to  an  internal  step  pressure. 


terms  of  O(s^),  whereas  the  DAA2  and  DAAi_2  relations,  (7.52)  and  (7.48),  include  terms  only  through 
O(s^). 

The  modal  solution  for  this  case  is  determined  by  using  tj  equal  to  po  in  equation  (7.36)  and  cti  equal 
to  2po.  No  (Ti  to  Ti  ratio  other  than  two  will  produce  a  stable  late-time  result.  Figure  8-5  presents  a 
comparison  of  radial  displacement  responses  for  the  exact  relation  (7.36),  the  DAA2  relation  (7.52), 
and  the  DAAi  relation  (7.46).  Figure  8-6  compares  the  meridional  (tangential)  displacement  solutions 
for  these  relations.  Agreement  between  the  exact  and  DAA2  modal  results  is  excellent. 

Finally,  Figures  8-7  and  8-8  compare  results  from  the  mixed-order  DAA,  DAAi_2,  (modal  and  BEM), 
with  exact  and  DAAi  •  We  see  that  the  essentially  coincident  modal  and  BEM  DAAi_2  histories  deviate 
modestly  from  the  exact  history. 

8.4  SPHERICAL  CAVITY/  HALF-SPACE/  STEP  PRESSURE  PROBLEM. 

The  third  problem  is  a  spherical  cavity  embedded  in  a  semi-infinite  elastic  medium  and  excited  by  an 
internal  step  pressure.  This  problem  is  very  difficult  for  doubly  asymptotic  approximations,  at  least 
as  formulated  to  date.  Only  the  late-time  approximations  have  terms  that  include  the  effect  of  the  free 
surface.  Specifically,  the  dynamic  free-surface  effect  is  manifested  only  in  the  half-space  terms  of  A, 
B,  and  f  [see  (5.7),  (5.8),  (5.12),  and  (5.14)].  The  reflection  of  the  waves  from  the  free  surface  occurs 
at  intermediate  times  if  the  cavity  is  fairly  close,  but  not  immediately  adjacent,  to  the  surface. 
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Tangential  Surface  Displacement  ^  Radial  Surface  Displacement 
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Exact,  DAAi,  and  DAA2  radial  displacements  of  a  spherical  cavity  in  an  infinite  elastic 
medium  subjected  to  n=  1  step  tractions. 


Figure  8-6.  Exact,  DAAi,  and,  DAA2  tangential  displacements  of  a  spherical  cavity  in  an  infinite 
elastic  medium  subjected  to  n  =  1  step  tractions. 
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Tangential  Surface  Displacement  ^  Radial  Surface  Displacement 
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Exact,  DAAi,  and  DAAi_2  radial  displacement  of  a  spherical  cavity  in  an  infinite  elastic 
medium  subjected  to  n  =  1  step  internal  tractions. 


Figure  8-8.  Exact,  DAAi,  and  DAAi_2  tangential  displacements  of  a  spherical  cavity  in  an  elastic 
medium  subjected  to  a  n  =  1  step  internal  tractions. 
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This  problem  does  not  possess  radial  symmetry,  and  does  not  posses  an  analytical  solution.  However, 
a  boundary-element  solution  based  on  numerical  inversion  of  Laplace  transforms  has  been  generated 
by  Manolis  and  Ahmad  (Manolis  and  Ahmad.,  1988),  and  an  analytical  solution  to  the  related  static 
problem  exists  (Bonefed,  1990).  In  terms  of  the  geometry  shown  in  Figure  8-9,  the  latter  solution  is 


=  Poa' 


R 


.5 


where 


r  r 


} 


Uz 


Poa^ 


2{l-2u) 


(-^  +  d) 


z  3^(2  +  d)^ 

P‘3  +  7?.5 


+  d)  {z  +  d)  1 

P'3  P*3  Jj 


P'  =  +  {z-  ZoY 

P*  =  +  {z  Zo)^. 


(8.7) 


Numerical  DAAi  and  numerical  inversion  solutions  for  this  problem  are  shown  in  Figures  8-10-8-12, 
along  with  the  late-time  static  asymptotes  given  by  (8.7);  the  physical  parameters  specified  are  the  same 
as  those  previously  given  for  the  infinite-domain  problem.  Figure  8-10  pertains  to  the  top  of  the  cavity, 
i.e.,  the  point  on  the  cavity  surface  closest  to  the  free  surface  of  the  elastic  half-space.  The  numerical 
inversion  solution  by  Manolis  and  Ahmad  shows  the  pressure  wave  reflected  from  the  free  surface,  at 
=  2.  The  effects  of  the  reflected  shear  and  Rayleigh  wave  at  =  2.73  and  —  3.4,  respec¬ 
tively,  complicate  the  response.  Figure  8-11  pertains  to  a  point  90°  around,  and  Figure  8-12  pertains  to 
the  bottom  of  the  cavity.  These  figures  show  that  the  DAAi  solutions  agree  with  the  numerical  inver¬ 
sion  solutions  at  early  time  and  appear  to  approach  the  correct  late-time  asymptotes.  Unfortunately, 
the  numerical  inversion  solutions  do  not  extend  far  enough  in  time  to  allow  a  completely  satisfactory 
comparison. 

We  also  generated  results  from  halfspace-D  A  A2  for  the  spherical  cavity  in  a  semi-infinite  elastic  domain 
subjected  to  an  internal  step  pressure.  We  show  the  DAA2  radial  displacements  at  three  locations 
on  the  spherical  surface,  6  =  0°,  ^  =  90°,  and  0  =  180°  in  Figures  8-13-8-15.  We  compare  these 
solutions  with  the  corresponding  numerical  inversion  solutions  (Manolis  and  Ahmad.,  1988)  and  the 
static  asymptotes  (t  =  00).  All  radial  displacements  have  been  divided  by  the  exact  static  solution  for 
a  spherical  cavity  in  an  infinite  elastic  domain,  u static  =  0- 

Radial  displacement  response  at  0  =  0°  (closest  to  free  surface)  shows  the  effect  of  the  dilatational 
wave  reflecting  from  the  free  surface  at  ^  =  2.0.  DAA2  radial  displacement  response  atO  =  90° 
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Figure  8-9.  Geometry  for  a  cavity  embedded  in  a  semi-infinite  elastic  medium. 


Figure  8-10.  Radial  displacement  response  of  a  step-pressurized  cavity 
in  a  semi-infinite  elastic  medium  {0  =  0°,  d  =  a). 
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Figure  8-11.  Radial  displacement  response  of  a  step-pressurized  cavity 
in  a  semi-infinite  elastic  medium  (9  =  90°,  d  =  2a). 
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Figure  8-12.  Radial  displacement  response  of  a  step-pressurized  cavity 
in  a  semi-infinite  elastic  medium  (9  =  180°,  d  =  3a). 
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Figure  8-13.  Radial  displacement  response  of  a  step-pressurized  cavity 

in  a  semi-infinite  elastic  medium  with  DAA2(0  =  0°,  cl  =  a). 


Figure  8-14.  Radial  displacement  response  of  a  step-pressurized  cavity 

in  a  semi-infinite  elastic  medium  with  DAA2  (6  =  90°, d  =  2a). 
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Figure  8-15.  Radial  displacement  response  of  a  step-pressurized  cavity  in  a  semi-infinite  elastic 
medium  with  DAA2  {9  =  180®,  d  =  3a). 


agrees  well  with  the  exact  solution  for  a  whole-space,  see  Section  8.2.  Radial  displacement  response 
at  ^  =  180°  exibits  the  greatest  overshoot  beyond  the  corresponding  late-time  asymptote. 

8.5  STEP  PRESSURE  ON  A  HALF-SPACE  WITH  WHOLE-SPACE  MATRICES. 

The  fourth  (two-dimensional)  canonical  problem  is  a  step  pressure  applied  to  a  portion  of  the  surface 
of  an  elastic  half-space  (Figure  8-16).  Because  of  the  spatial  discontinuity  of  the  load,  a  relief  wave 
travels  from  each  discontinuity  toward  r  =  0.  The  problem  has  been  studied  analytically  by 
Eason  (Eason,  1966).  The  surface  is  flat  so  there  is  no  curvature,  a  =  1.61,  and  the  material  properties 
are  the  same  as  those  used  previously.  The  surface  is  discretized  using  36  elements.  We  terminated  the 
model  mesh  at  two  different  radii,  ^  =  1.0  (Manoiis  and  Beskos,  1988)  and  —  =  16.0.  We  compare 
results  for  the  normal  displacement  at  r  =  0  for  the  two  boundary  radii  in  Figure  8-17.  Eason’s  analyt¬ 
ical  result  is  suspect  at  ^  =  1.8;  the  discontinuous  displacement  is  probably  caused  by  the  inversion 
method  used  (Kim  and  Soedel,  1988).  Note  that  the  result  from  the  second-order  whole-space  DAA 
does  have  a  smooth  response.  The  smaller  radius  mesh  begins  to  diverge  significantly  at  ^  =  5.0. 
The  second-order  DAA  response  is  linear  until  ^  =  1.0,  when  a  dilatational  wave  reaches  the  center 
(r  =  0)  from  the  edge  of  the  loaded  area.  When  ^  >  1.0  the  second-order  DAA  strongly  resembles 
a  first-order  DAA,  i.e.  it  produces  no  elastic  overshoot.  This  is  because  the  surface  is  flat  so  the  cur¬ 
vature  terms  in  the  second-order  early-time  approximations  are  zero  for  this  problem.  As  expected  the 
late-time  (static)  result  is  accurately  found. 
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f/Pa 


Figure  8-18.  Normal  displacement  at  r  =  0  for  normal  step-load  applied  to  an  elastic  half-space  with 
half-space  DAA2. 


8.6  STEP  PRESSURE  ON  A  HALF-SPACE  WITH  HALF-SPACE  MATRICES. 

The  last  comparison  also  pertains  to  the  half-space  problem.  The  discretatization  for  this  problem  con¬ 
sisted  of  36  elements.  Figure  8-18  shows  normal-displacement  response  at  r  =  0  from  the  half-space 
DAA2  along  with  Eason’s  analytical  result.  The  half-space  DAA2  response  is  almost  identical  to  the 
whole-space  DAA2  response  discussed  previously,  but  was  computed  in  about  one-fourth  the  time  be¬ 
cause  no  boundary  elements  were  needed  for  the  unloaded  surface  area. 
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SECTION  9 


CONCLUSION 


This  report  documents  a  systematic  formulation  and  an  analytical  and  numerical  evaluation  of  first- 
order  and  second-order  doubly  asymptotic  approximations  for  computational  boundaries  in  the  solution 
of  transient  soil-structure  interaction  problems.  Specifically, 

1.  The  early-time  (high-frequency)  approximations  for  an  acoustic  medium,  systematically  formu¬ 
lated  by  Felippa  (Felippa,  1980b),  were  extended  to  isotropic  elastic  domains. 

2.  Ray  elastodynamics,  developed  by  H.  B.  Keller  and  J.  B.  Keller  (Keller,  1958),  (Keller,  1964), 
(Ahluwalia  et  aL,  1969),  was  also  employed  in  the  formulation  of  the  early-time  approximations. 

3.  The  first-  and  second-order  late-time  (low-frequency)  approximations  for  infinite  and  semi-infinite 
elastic  media  were  systematically  formulated  from  the  elastodynamic  boundary-integral  equa¬ 
tions  of  Cruse  and  Rizzo  (Cruse  and  Rizzo,  1968)  and  the  half-space  operators  of  Banerjee  and 
Mamoon  (Banerjee  and  Mamoon,  1990). 

4.  The  first-order  DAA’s  for  both  infinite  and  semi-infinite  isotropic  elastic  media  were  systemati¬ 
cally  formulated;  the  operator  matching  method  of  Nicholas- Vullierme  (Nicholas- Vullierme,  1991) 
and  of  Geers  and  Zhang  (Geers  and  Zhang,  1991)  was  used  to  derive  these  DAA’s,  which  support 
the  heuristic  formulation  of  Underwood  and  Geers  (Underwood  and  Geers,  1981). 

5.  The  second-order  DAA’s  for  both  infinite  and  semi-infinite  isotropic  elastic  media  were  system¬ 
atically  formulated  using  the  operator-matching  method. 

6.  Mixed-order  DAA’s  were  formulated  for  both  inifinite  and  semi-infinite  isotropic  elastic  media; 
each  combines  a  first-order  low-frequency  approximation  with  a  second-order  high-frequency 
approximation. 

7 .  Finite-element  discretization  was  introduced  to  configure  the  D AA  formulas  for  boundary-element 
solution;  pertinent  computer  programs  were  written  to  generate  numerical  results. 

8.  Exact,  DAAi,  DAA2  and  DAAi_2  modal  computational  boundary  relations  were  derived  for  di- 
latational,  rotational  and  translational  motions  of  a  spherical  boundary  in  an  infinite  isotropic 
elastic  medium. 

9.  Boundary-element  DAAi  and  DAAi  _2  results  for  suddenly  pressurized  spherical  cavities  embed¬ 
ded  in  infinite  and  semi-infinite  elastic  media  were  generated;  the  DAA  results  were  compared 
with  corresponding  analytical  results  from  Timoshenko  and  Goodier 

(Timoshenko  and  Goodier,  1970),  Manolis  and  Ahmad  (Manolis  and  Ahmad.,  1988), 

Bonefed  (Bonefed,  1990),  and  with  analytical  results  generated  herin. 


66 


10.  Boundary-element  DAA2  results  for  a  pressure  load  applied  suddenly  to  a  portion  of  the  surface 
of  an  elastic  semi-infinite  domain  were  generated;  these  were  compared  with  results  reported  by 
Eason  (Eason,  1966). 

The  principal  conclusions  reached  in  this  study  are: 

1 .  First-order  doubly  asymptotic  approximations  are  useful  as  computational  boundaries  for  elastic 
media,  but  they  cannot  predict  the  response  overshoot  that  characterizes  many  problems. 

2.  Second-order  doubly  asymptotic  approximations  are  quite  useful  as  computational  boundaries 
for  elastic  media,  successfully  predicting  response  overshoot. 

3.  Half-space  DAA’s  are  much  more  efficient  than  whole-space  DAA’s  for  half-space  problems  be¬ 
cause  the  former  do  not  require  discretization  of  unloaded  free-surface  areas. 

4.  The  whole-space  DAA2  is  exact  for  dilatational  and  rotational  motion  of  a  spherical  surface;  for 
the  translational  cases  studied,  the  modal  DAA2  results  agree  perfectly  with  the  corresponding 
exact  results. 

5.  The  mixed-order  whole-space  DAA  (DAAi_2)  is  exact  for  dilatational  and  rotational  motion  of 
a  spherical  surface;  it  is  less  accurate  than  the  whole-space  DAA2  for  translational  motion. 
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APPENDIX 


TENSOR  OPERATORS  FOR  HALF  SPACE 


The  following  are  the  first  group  of  U  operators  that  pertain  to  the  cylindrical  coordinate  system. 
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